IISc-CTS-12/99 
ZU-TH 10/00 
BUTP-99/33 



Roy equation analysis of nn scattering 

May 30, 2000 

B. Ananthanarayan a , G. Colangelo 6 , J. Gasser c and H. Leutwyler c 



a Centre for Theoretical Studies, Indian Institute of Science 
Bangalore, 560 012 India 
b Institute for Theoretical Physics, University of Zurich 
Winterthurerstr. 190, CH-8057 Zurich, Switzerland 
c Institute for Theoretical Physics, University of Bern 
Sidlerstr. 5, CH-3012 Bern, Switzerland 



Abstract 

We analyze the Roy equations for the lowest partial waves of elastic 
7T7T scattering. In the first part of the paper, we review the mathematical 
properties of these equations as well as their phenomenological applica- 
tions. In particular, the experimental situation concerning the contribu- 
tions from intermediate energies and the evaluation of the driving terms 
are discussed in detail. We then demonstrate that the two S'-wave scatter- 
ing lengths clq and ag are the essential parameters in the low energy region: 
Once these are known, the available experimental information determines 
the behaviour near threshold to within remarkably small uncertainties. An 
explicit numerical representation for the energy dependence of the S- and 
P-waves is given and it is shown that the threshold parameters of the D- 
and F- waves are also fixed very sharply in terms of aj] and Oq. In agree- 
ment with earlier work, which is reviewed in some detail, we find that the 
Roy equations admit physically acceptable solutions only within a band of 
the (g,q,(Iq) plane. We show that the data on the reactions e + e~ — » tttt and 
r —* 7r 7r v reduce the width of this band quite significantly. Furthermore, 
we discuss the relevance of the decay K — > tt it e v in restricting the allowed 
range of a{j , preparing the grounds for an analysis of the forthcoming pre- 
cision data on this decay and on pionic atoms. We expect these to reduce 
the uncertainties in the two basic low energy parameters very substantially, 
so that a meaningful test of the chiral perturbation theory predictions will 
become possible. 
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1 Introduction 

The present paper deals with the properties of the tttt scattering amplitude in 
the low energy region. Our analysis relies on a set of dispersion relations for the 
partial wave amplitudes due to Roy [|TJ . These equations involve two subtraction 
constants, which may be identified with the S'-wave scattering lengths, a[j and Oq. 
We demonstrate that the subtraction constants represent the essential parameters 
in the low energy region - once these are known, the Roy equations allow us 
to calculate the partial waves in terms of the available data, to within small 
uncertainties. Given the strong dominance of the two S'-waves and of the P-wave, 
it makes sense to solve the equations only for these, using experimental as well as 
theoretical information to determine the contributions from higher energies and 
from the higher partial waves. More specifically, we solve the relevant integral 
equations on the interval 2M n < yfs < 0.8 GeV. One of the main results of this 
work is an accurate numerical representation of the S- and P-waves for a given 
pair of scattering lengths a[j and Oq. 

Before describing the outline of the present paper, we review previous work 
concerning the Roy equations. Roy's representation |IJ for the partial wave am- 
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plitudes tj of elastic tttc scattering reads 



2 00 poo 

ti( s ) = kj( s ) + YY / ds'Kfi, ( s ,s')im4(s') , (1.1; 

/4M2 



I'=0 £'=0 



where I and i denote isospin and angular momentum, respectively and kj(s) is 
the partial wave projection of the subtraction term. It shows up only in the S- 
and P- waves, 

His) = 4 5° + S -^f (2a» - 5al) (i 5* 6° + i S{ 5] - \ 5{ <5°) . (1.2) 

The kernels K^j, (s, s') are explicitly known functions (see appendix A). They 
contain a diagonal, singular Cauchy kernel that generates the right hand cut 
in the partial wave amplitudes, as well as a logarithmically singular piece that 
accounts for the left hand cut. The validity of these equations has rigorously 
been established on the interval — 4M^ < s < 60M^. 

The relations ( |1.1| ) are consequences of the analyticity properties of the irn 
scattering amplitude, of the Froissart bound and of crossing symmetry. Com- 
bined with unitarity, the Roy equations amount to an infinite system of coupled, 
singular integral equations for the phase shifts. The integration is split into a 
low energy interval 4M^ < s' < Sq and a remainder, sq < s' < 00. We refer to 
So as the matching point, which is chosen somewhere in the range where the Roy 
equations are valid. The two S'-wave scattering lengths, the elasticity parameters 
below the matching point and the imaginary parts above that point are treated 
as an externally assigned input. The mathematical problem consists in solving 
Roy's integral equations with this input. 

Soon after the original article of Roy had appeared, extensive phenomeno- 
logical applications were performed @-||, resulting in a detailed analysis and 
exploitation of the then available experimental data on tot scattering. For a re- 
cent review of those results, we refer the reader to the article by Morgan and 
Pennington 0. Parallel to these phenomenological applications, the very struc- 
ture of the Roy equations was investigated. In O], a family of partial wave 
equations was derived, on the basis of manifestly crossing symmetric dispersion 
relations in the variables st + tu + u s and stu. Each set in this family is valid 
in an interval Sq < s < si, and the union of these intervals covers the domain 
— 28M^ < Re s < 125. 3M^ (for a recent application of these dispersion relations, 
see ||12|| ). Using hyperbolae in the plane of the above variables, Auberson and 
Epele [0 proved the existence of partial wave equations up to Re s = 165M^. 



Furthermore, the manifold of solutions of Roy's equations was investigated, in 



the single channel [14|— [16| as well as in the coupled channel case [17|. In the 
late seventies, Pool |T8| provided a proof that the original, infinite set of integral 
equations does have at least one solution for < 4.8 M^, provided that the 



driving terms are not too large, see also [19]. Heemskerk and Pool also examined 
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numerically the solutions of the Roy equations, both by solving the N equation 
[|l9|j and by using an iterative method |20| . 

It emerged from these investigations that - for a given input of .S-wave scat- 
tering lengths, elasticity parameters and imaginary parts - there are in general 
many possible solutions to the Roy equations. This non-uniqueness is due to the 



singular Cauchy kernel on the right hand side of (|1 . 1| ) . In order to investigate the 



uniqueness properties of the Roy system, one may - in a first step - keep only 
this part of the kernels, so that the integral equations decouple: one is left with a 
single channel problem, that is a single partial wave, which, moreover, does not 
have a left hand cut. This mathematical problem was examined by Pomponiu 
and Wanders, who also studied the effects due to the presence of a left hand cut 
H|j . Investigating the infinitesimal neighbourhood of a given solution, they found 



that the multiplicity of the solution increases by one whenever the value of the 
phase shift at the matching point goes through a multiple of 7r/2. Note that the 
situation for the usual partial wave equation is different: There, the number of 
parameters in general increases by two whenever the phase shift at infinity passes 



through a positive integer multiple of ir, see for instance |21|, Q an d references 
cited therein. 

After 1980, interest in the Roy equations waned, until recently. For instance, 
in refs. [|3| these equations are used to analyze the threshold parameters for 
the higher partial waves, relying on the approach of Basdevant, Froggatt and 
Petersen [[5], ||. The uncertainties in the values of a[] and cZq are reexamined in 
refs. pi[ . In recent years, it has become increasingly clear, however, that a new 
analysis of the itit scattering amplitude at low energies is urgently needed. New 
K e4: experiments and a measurement of the combination a[j — Oq based on the 



decay of pionic atoms are under way p5f-p9f. It is expected that these will 



significantly reduce the uncertainties inherent in the data underlying previous 
Roy equation studies, provided the structure of these equations can be brought 
under firm control. For this reason, the one-channel problem has been revisited 
in great detail in a recent publication [30], while the role of the input in Roy's 
equations is discussed in ref. |HJ . 

The main reason for performing an improved determination of the irn scat- 
tering amplitude is that this will allow us to test one of the basic properties of 
QCD, namely the occurrence of an approximate, spontaneously broken symme- 
try: The symmetry leads to a sharp prediction for the two S'-wave scattering 



lengths ||32||— 1|40|| . The prediction relies on the standard hypothesis, according to 
which the quark condensate is the leading order parameter of the spontaneously 
broken symmetry. Hence an accurate test of the prediction would allow us to ver- 
ify or falsify that hypothesis ||34|| . First steps in this program have already been 
performed |5S[|-[p9[l. However, in the present paper, we do not discuss this issue. 
We follow the phenomenological path and ignore the constraints imposed by chi- 
ral symmetry altogether, in order not to bias the data analysis with theoretical 
prejudice. In a future publication, we intend to match the chiral perturbation 
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theory representation of the scattering amplitude to two loops ||0 with the phe- 
nomenological one obtained in the present work. 

Finally, we describe the content of the present paper. Our notation is specified 
in section 2. Sections 3 and 4 contain a discussion of the background amplitude 
and of the driving terms, which account for the contributions from the higher 
partial waves and from the high-energy region. As is recalled in section 5, uni- 
tarity leads to a set of three singular integral equations for the two S'-waves and 
for the P-wave. The uniqueness properties of the solutions to these equations are 
discussed in section 6, while section 7 contains a description of the experimental 
input used for energies between 0.8 and 2 GeV. In particular we also discuss 
the information concerning the P-wave phase shift, obtained on the basis of the 
e + e~ — > 7T7T and r — > tt-kv data. In section 8, we describe the method used to 
solve the integral equations for a given input. The resulting universal band in the 
(a[],aQ) plane is discussed in section 9, where we show that, in the region below 
0.8 GeV, any point in this band leads to a decent numerical solution for the three 
lowest partial waves. As discussed in section 10, however, the behaviour of the 
solutions above that energy is consistent with the input used for the imaginary 
parts only in part of the universal band - approximately the same region of the 
(clq,clI) plane, where the Olsson sum rule is obeyed (section 11). The solutions 
are compared with available experimental data in section 12, and in section 13, 
we draw our conclusions concerning the allowed range of a[j and a^. The other 
threshold parameters can be determined quite accurately in terms of these two. 
The outcome of our numerical evaluation of the scattering lengths and effective 
ranges of the lowest six partial waves as functions of a[j and Oq is given in section 
14, while in section 15, we describe our results for the values of the phase shifts 
relevant for K — > tttt. Section 16 contains a comparison with earlier work. A 
summary and concluding remarks are given in section 17. 

In appendix A we describe some properties of the Roy kernels, which are 
extensively used in our work. The background from the higher partial waves 
and from the high energy tail of the dispersion integrals is discussed in detail 
in appendix B. In particular, we show that the constraints imposed by crossing 
symmetry reduce the uncertainties in the background, so that the driving terms 
can be evaluated in a reliable manner. In appendix C we discuss sum rules 
connected with the asymptotic behaviour of the amplitude and show that these 
relate the imaginary part of the P-wave to the one of the higher partial waves, 
thereby offering a sensitive test of our framework. Explicit numerical solutions of 
the Roy equations are given in appendix D and, in appendix E, we recall the main 
features of the well-known Lovelace-Shapiro- Veneziano model, which provides a 
useful guide for the analysis of the asymptotic contributions. 
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2 Scattering amplitude 

We consider elastic Tin scattering in the framework of QCD and restrict our 
analysis to the isospin symmetry limit, where the masses of the up and down 
quarks are taken equal and the e.m. interaction is ignored^. In this case, the 
scattering process is described by a single Lorentz invariant amplitude A(s,t,u), 

(Tr d (pA)TT c (p 3 ) out|7T a (p 1 )7r 6 (p 2 ) in) = 5 fi + 

(2n)H 6 4 {P f - P l ){5 ab 5 cd A(s, t, u) + 5 ac 5 bd A(t, u, s) + 5 ad 5 bc A(u, s,t)} . 

The amplitude only depends on the Mandelstam variables s, t, u, which are 
constrained by s + t + u = AM 2 . Moreover, crossing symmetry implies 

A(s, t, u) = A(s, u, t) . 

The s-channel isospin components of the amplitude are given by 

T°(s, t) = 3A(s, t, u) + A(t, u, s) + A{u, s, t) , 
T 1 (s,t)=A(t,u,s) - A{u,s,t) , (2.1) 
T 2 (s,t) =A(t,u,s) +A(u,s,t) . 

In our normalization, the partial wave decomposition reads 



T 7 (s,t)=32 7r^(2£+ l)P e (1 + 



2* \ T, s 



2ia(s 



s - 4MJ / " 

±- «»«<•) -1} , (2.2) 



/ 4M 2 

The threshold parameters are the coefficients of the expansion 

Ret I e {s) = q 2e {a I e +q 2 b I e +q 4 c I e + ...} , (2.3) 

with s = A(Ml + q 2 ). 

The isospin amplitudes T = (T°,T l ,T 2 ) obey fixed-t dispersion relations, 
valid in the interval -28M 2 < t < AM 2 f4l| . As shown by Roy 0, these can be 
written in the formQ 

f(s, t) = {AMI)- 1 (sl + tC st + u C su ) f(AM 2 n , 0) (2.4) 

+ / ds' g 2 (s,t, s')lmT(s' ,0) + / ds' g s (s,t, s')lm.T(s' ,t) . 
JiMl Jam$ 



our numerical work, we identify the value of M v with the mass of the charged pion. 
2 For an explicit representation of the kernels 172(5, i, s'), <?3(s, t, s') and of the crossing ma- 
trices C s t, C su , we refer to appendix A. 
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The subtraction term is fixed by the S'-wave scattering lengths: 

f(AM 2 ,0)=32n(a° ,0,a 2 ). 

The Roy equations ( |1.1| ) represent the partial wave projections of eq. (|2. 1|). 
Since the partial wave expansion of the absorptive parts converges in the large 
Lehmann-Martin ellipse, these equations are rigorously valid in the interval 
— AM 2 < s < 60M 2 . If the scattering amplitude obeys Mandelstam analyti- 
city, the fixed-t dispersion relations can be shown to hold for — 32M 2 < t < AM 2 
and the Roy equations are then also valid in a larger domain: —AM 2 < s < 68M 2 
(for a review, see |32|). In fact, as we mentioned in the introduction, the range 
of validity can be extended even further [[LI], JT^], so that Roy equations could 
be used to study the behaviour of the partial waves above a/68 M w = 1.15 GeV, 
where the uncertainties in the data are still considerable. In the following, how- 
ever, we focus on the low energy region. We assume Mandelstam analyticity and 
analyze the Roy equations in the interval from threshold to 

si = 68M 2 , v^7= 1.15 GeV . 

3 Background amplitude 

The dispersion relation (|2.4Q shows that, at low energies, the scattering amplitude 
is fully determined by the imaginary parts of the partial waves in the physical 
region, except for the two subtraction constants a®, a 2 ,. In view of the two sub- 
tractions, the dispersion integrals converge rapidly. In the region between 0.8 and 
2 GeV, the available phase shift analyses provide a rather detailed description of 
the imaginary parts of the various partial waves. Our analysis of the Roy equa- 
tions allows us to extend this description down to threshold. For small values of s 
and t, the contributions to the dispersion integrals from the region above 2 GeV 
are very small. We will rely on Regge asymptotics to estimate these. In the fol- 
lowing, we split the interval of integration into a low energy part (AM 2 < s' < S2) 
and a high energy tail (S2 < s' < 00), with 

= 2 GeV , s 2 = 205.3 M 2 . 

For small values of s and t, the scattering amplitude T(s, t) is dominated by 
the contributions from the subtraction constants and from the low energy part 
of the dispersion integral over the imaginary parts of the S- and P-waves. We 
denote this part of the amplitude by T(s,t) SP . The corresponding contribution 
to the partial waves is given by 

2 1 

= + E E [\ds'K I u i :(s,s')lmt i ;(s') . (3.1) 
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The remainder of the partial wave amplitude, 

2 oo 

d\{s) =EE f\ ds'K»'(s, s') lm4(s') (3.2) 

2 °° POO 

i'=oe'=o Js2 

is called the driving term. It accounts for those contributions to the r.h.s. of the 
Roy equations that arise from the imaginary parts of the waves with £ = 2,3,... 
and in addition also contains those generated by the imaginary parts of the S- 
and P-waves above 2 GeV. By construction, we have 

t\{s)=t\{s) SP + d I t {s) . (3.3) 

For the scattering amplitude, the corresponding decomposition reads 

f(s,t)=f(s,t) SP + f(s,t) d . (3.4) 

We refer to T(s, t)d as the background amplitude. 

The contribution from the imaginary parts of the S- and P-waves turns out 
to be crossing symmetric by itself. In this sense, crossing symmetry does not 
constrain the imaginary parts of the S- and P-wavesQ. The symmetry can be 
exhibited explicitly by representing the three components of the vector T(s, t) SP as 
the isospin projections of a single amplitude A(s, t, u) sp that is even with respect 
to the exchange of t and u. The explicit expression involves three functions of a 
single variable [IT, 3U|: 

A(s, t, u) SP = 32tt [lW°{s) + f (s - u)W\t) + |(s - t)^ 1 ^) 

+ \W 2 (t) + \W 2 (u)-\W 2 (s)} . (3.5) 

These are determined by the imaginary parts of the S- and P-waves and by the 
two subtraction constants a[}, a^: 

,,, , a° s s(s-AM 2 ) y«a ds 1 Imtg(s') 



1 J AMI 7T Ja 



4M2 s'(s' -4M%)(s' - s) 



„,-,, , s f s 2 ds' ImtUs') 



w2(s a 2 s , s(s-4M3 r a tfe'Im^(s') 



4M2 7T J 4 Af2 s'(s' - 4Af2)(s' - s) 

The representation 

A(s, t, u) = A(s, t, u) SP + A(s, t, u) d (3.7) 

3 The asymptotic behaviour of the scattering amplitude does tie the imaginary part of the 
P-wave to the contributions from the higher partial waves, see appendix C.l. 
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yields a manifestly crossing symmetric decomposition of the scattering amplitude 
into a leading term generated by the imaginary parts of the S- and P-waves at 
energies below S2 and a background, arising from the imaginary parts of the 
higher partial waves and from the high energy tail of the dispersion integrals. 

4 Driving terms 

In the present paper, we restrict ourselves to an analysis of the Roy equations 
for the S- and P- waves, which dominate the behaviour at low energies. The 
background amplitude only generates small corrections, which can be worked 
out on the basis of the available experimental information. The calculation is 
described in detail in appendix [B]. In particular, we show that crossing symmetry 
implies a strong constraint on the asymptotic contributions. 

The resulting numerical values for the driving terms are well described by 
polynomials in s, or, equivalently, in the square of the center of mass momentum 
q 2 = ~(s — 4M 2 ). By definition, the driving terms vanish at threshold, so that 
the polynomials do not contain g-independent terms. In view of their relevance 
in the evaluation of the threshold parameters, we fix the coefficients of the terms 
proportional to q 2 with the derivatives at threshold and also pin down the term 
of order g 4 in the P-wave, such that it correctly accounts for the background 
contribution to the effective range of this partial wave. The remaining coefficients 
of the polynomial are obtained from a fit on the interval from threshold to s±. 
The explicit result reads 

d° (s) = 0.116 q 2 + 4.79 g 4 - 4.09 g 6 + 2.69 g 8 , 

d\(s) = 0.00021 q 2 + 0.038 g 4 + 0.94 g 6 - 1.21 q 8 , (4.1) 

d 2 (s) = 0.0447 q 2 + 1.59 g 4 - 6.26 g 6 + 5.94 g 8 , 

where g is taken in GeV units (the range 4M 2 < s < 68M 2 corresponds to 
< q < 0.56 GeV). The driving term of the 1 = S-wave is larger than the 
other two by an order of magnitude. It is dominated almost entirely by the 
contribution from the P-wave with 1 = 0. In d\(s), the P- and P-waves nearly 
cancel, so that the main contributions arise from the region above 2 GeV. The 
term rfg(s) picks up small contributions both from low energies and from the 
asymptotic domain. The above polynomials are shown as full lines in fig. [TJ. The 
shaded regions represent the uncertainties of the result, which may be represented 
as d{(s) ± e|(s), with 

eg(s) = 0.008 q 2 + 0.31 g 4 - 0.33 g 6 + 0.41 g 8 , 

e^(s)=0.002g 2 + 0.06g 4 -0.17g 6 + 0.21g 8 , (4.2) 

e 2 (s) = 0.005 g 2 + 0.20 g 4 - 0.32 g 6 + 0.39 g 8 . 

Above threshold, the error bars in d|j(s), d\(s) and dl(s) roughly correspond to 
6%, 1% and 4% of cZq(s) , respectively. 
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Figure 1: Driving terms versus energy in GeV. The full lines show the re- 
sult of the calculation described in appendix B. The shaded regions indicate the 
uncertainties associated with the input of that calculation. The dashed curves 
represent the contributions from the D- and F-waves below 2 GeV. 
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As far as %(s) is concerned, our result roughly agrees with earlier calculations 
[[|, ||. Our values for d\(s) and g?q(s), however, are much smaller. The bulk of 
the difference is of purely kinematic origin: The values taken for S2 are different. 
While we are working with yJ~S2 = 2 GeV, the values used in refs. 0] and || are 
\^53M n ~ 1 GeV and VHO M n ~ 1.5 GeV, respectively. The value of S2 enters 
the definition of the driving terms in eq. (3.2) as the lower limit of the integration 
over the imaginary parts of the S- and P- waves. We have checked that, once this 
difference in the range of integration is accounted for, the driving terms given 
in these references are consistent with the above representation. Note however, 
that our uncertainties are considerably smaller, and we do rely on this accuracy 
in the following. It then matters that not only the range of integration, but also 
the integrands used in |3], Q differ from ours: In these references, it is assumed 
that, above the value taken for S2, the behaviour of the S- and P-wave imaginary 
parts is adequately described by a Regge representation. 

The difference between such a picture and our representation for the back- 
ground amplitude is best illustrated with the simple model used in the early 
literature, where the asymptotic region is described by a Pomeron term with 
a tot — 20 mb and a contribution from the p -/-trajectory, taken from the Lovelace- 
Shapiro- Veneziano model (appendix [E]). As discussed in detail in appendix |B.4 



the assumption that an asymptotic behaviour of this type sets in early is in con- 



flict with crossing symmetry [ 4"3"| . In particular, the model overestimates the 
contribution to the driving terms from the region above 1.5 GeV, roughly by a 
factor of two. Either the value of a tot or the residue of the leading Regge trajec- 
tory or both must be reduced in order for the model not to violate the sum rule 
( |B.(j| ). The manner in which the asymptotic contribution is split into one from 
the Pomeron and one from the leading Regge trajectory is not crucial. For any 
reasonable partition that obeys the sum rule (|B.6|) , the outcome for the driving 
terms is approximately the same. The result for d\(s) and cIq(s) is considerably 
smaller than what is expected from the above model. The leading term cIq(s), 
on the other hand, is dominated by the resonance ^2(1275) and is therefore not 
sensitive to the behaviour of the imaginary parts in the region above 1.5 GeV. 



5 Roy equations as integral equations 

Once the driving terms are pinned down, the Roy equations for the S- and P- 
waves express the real parts of the partial waves in terms of the S'-wave scattering 
lengths and of a principal value integral over their imaginary parts from 4M^ to 
S2- Unitarity implies that, in the elastic domain 4M^ < s < 16M%, the real 
and imaginary parts of the partial wave amplitudes are determined by a single 
real parameter, the phase shift. If we were to restrict ourselves to the elastic 
region, setting s 2 = 16M^, the Roy equations would amount to a set of coupled, 
nonlinear singular integral equations for the phase shifts. We may extend this 
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range, provided the elasticity parameters 77/(5) are known. On the other hand, 
since the Roy equations do not constrain the behaviour of the partial waves 
for s > 68 M%, the integrals occurring on the r.h.s. of these equations can be 
evaluated only if the imaginary parts in that region are known, together with 
the subtraction constants Oq, Oq, which also represent parameters to be assigned 
externally. 

In the present paper, we do not solve the Roy equations in their full domain 
of validity, but use a smaller interval, 4M^ < s < s . The reason why it is 
advantageous to use a value of so below the mathematical upper limit, so < si, 
is that the Roy equations in general admit more than one solution. As will be 
discussed in detail in section || the solution does become unique if the value of 
So is chosen between the p mass and the energy where the 1 = S-w&ve phase 
passes through tt/2 - this happens around 0.86 GeV. In the following, we use 

v^o" = 0.8 GeV , s = 32.9 Ml . 

In the variable s, our matching point is nearly at the center of the interval between 
threshold and s\ = 68 M%. We are thus solving the Roy equations on the lower 
half of their range of validity, using the upper half to check the consistency of 
the solutions so obtained (section |n]). Our results are not sensitive to the precise 
value taken for s (section |9|). 

The Roy equations for the S- and P-waves may be rewritten in the form 

Re^(s) = kj{s)+ + ds'K^(s, s') Imt°(s')+ f ds'K^(s, s') lmt\(s') 

J4M% J4M% 

+ -f S \ ds'K^(s, s') Imtl(s') + fl(s) + d{(a) , (5.1) 

J AMI 

where I and £ take only the values (I,£) =(0,0), (1,1) and (2,0). The bar across 
the integral sign denotes the principal value integral. The functions //(s) contain 
the part of the dispersive integrals over the three lowest partial waves that comes 
from the region between so and S2, where we are using experimental data as 
input. They are defined as 

ft is) = E E F ds>K»'(s } s') lines') . (5.2) 

The experimental input used to evaluate these integrals will be discussed in sec- 
tion |7|, together with the one for the elasticity parameters of the S- and P-waves. 

One of the main tasks we are faced with is the construction of the numerical 
solution of the integral equations (|5.1|) in the interval 4M^ < s < sq, for a given 
input {a[], Oq, //, rjj, dj}. Once a solution is known, the real part of the amplitude 
can be calculated with these equations, also in the region Sq < s < si. 
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6 On the uniqueness of the solution 



The literature concerning the mathematical structure of the Roy equations was 
reviewed in the introduction. In the following, we first discuss the situation for 
the single channel case - which is simpler, but clearly shows the salient features - 
and then describe the generalization to the three channel problem we are actually 
faced with. For a detailed analysis, we refer the reader to two recent papers on 
the subject |3(| |3l| and the references quoted therein. 



6.1 Roy's integral equation in the one-channel case 



If we keep only the diagonal, singular Cauchy kernel in (LI), the partial wave 
relations decouple, and the left hand cut in the amplitudes disappears. Each one 
of the three partial wave amplitudes then obeys the following conditions: 

i) In the interval between the threshold s = 4M% and the matching point s = sq, 
the real part is given by a dispersion relation 

, 1 r 00 , Imt(s') 
Ret(s) = a + (s-4M*)-l ds' — ,■,,;/, r- 6.1 

W nJ IT J4M* (s' - AM*) ( S ' - s) 

ii) Above Sq, the imaginary part lmt(s) is a given input function 

]mt(s) = A(s), s>s Q . (6.2) 

iii) For simplicity, we take the matching point in the elastic region, so that 

t( s ) = -^-e i5(s) sin S(s) , AMI < s < s , (6.3) 
a(s) 

where 5(s) is real and vanishes at threshold. We refer the reader to [3(J for a 



precise formulation of the regularity properties required from the amplitude and 
from the input absorptive part. As a minimal condition, we must require 

lim Im t(s) = A(s ) . (6.4) 

s/s 

Otherwise, the principal value integral does not exist at the matching point. 

Equations (|6.1| )- (|6.4j ) constitute the mathematical problem we are faced with 
in this case: Determine the amplitudes t(s) that verify these equations for a given 
input of scattering length a and absorptive part A(s). Once a solution is known, 
the real part of the amplitude above so is obtained from the dispersion relation 
( |6.1|) , and t(s) is then defined on 4M^ < s < oo. The following points summarize 
the results relevant in our context: 

1. Elastic unitarity reduces the problem to the determination of the real func- 
tion 5(s), defined in the interval 4M^ < s < Sq. The amplitude t(s) is then 
obtained from 
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a)m = b) m = 1 c)m = 2 

Figure 2: Boundary conditions on the phase S(so) for solving Roy's integral 
equation. Figs. a,b,c represent the cases < 5(sq) < ir/2, tt/2 < 8(sq) < n and 
7r < 8(sq) < 3n/2, respectively. In fig. c, the phase winds around the Argand 
circle slightly more than once. 



2. A given input {a, A(s)} does not, in general, fix the solution uniquely - in 
addition, the value of the phase at the matching point plays an important 
role. Indeed, let t(s) be a solution and suppose first that the phase at 
the matching point is positive. For < 8(s ) < tc/2, the infinitesimal 
neighbourhood of t(s) does not contain further solutions. For 8(sq) > n/2, 
however, the neighbourhood contains an m-parameter family of solutions. 
The integer m is determined by the value of the phase at the matching 
point ([x] is the largest integer not exceeding x): 

~25(s Q ) 
m = 

7T 

For a monotonically increasing phase, the index m counts the number of 
times 8(s) goes through multiples of 7r/2 as s varies from threshold to the 
matching point. We illustrate the situation for m = 0, 1, 2 in figure 0. 

3. If the value of the phase at the matching point is negative, the problem 
does not in general have a solution. In order for the problem to be soluble 
at all, the input must be tuned. For — n/2 < 8(s ) < 0, for instance, we 
may keep the absorptive part A(s) as it is, but tune the scattering length a. 
This situation may be characterized by m = —1: Instead of having a family 
of solutions containing free parameters, the input is subject to a constraint. 
Once a solution does exist, it is unique in the sense that the infinitesimal 
neighbourhood does not contain further solutions. 

4. Consider now the case displayed in fig. ^a, where the phase at the matching 
point is below tt/2. This corresponds to the situation encountered in the 
coupled channel case, for our choice of the matching point. According to the 



(6.5) 
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above statements, a given input {a,A(s)} then generates a locally unique 
solution - if a solution exists at all. We take it that uniqueness also holds 
globally, see ]T3 . 



The solution may be constructed in the following manner: Consider a family of 
unitary amplitudes, parametrized through Ci, . . . ,c n . For any given amplitude, 
evaluate the right and left hand sides of eq. (|6.1| ) and calculate the square of the 
difference at N points in the interval 4M^ < s < sq. Finally, minimize the sum 
of these squares by choosing ci, . . . , c n accordingly. Since the solution is unique, 
it suffices to find one with this method - it is then the only one. 

6.2 Cusps 

In general, the solutions are not regular at the matching point, but have a cusp 
(branch point) there: 5(s) = 5(s Q )+C(s — s) 7 + . . . , with 7 > 0. The phenomenon 
arises from our formulation of the problem - the physical amplitude is regular 
there. We conclude that, even if a mathematical solution can be constructed for 
a given input {a,A(s)}, it will in general not be acceptable physically, because 
it contains a fictitious singularity at the matching point. The behaviour of the 
phase is sensitive to the value of the exponent: If 7 is close to 1, the discontinuity 
in the derivative is barely visible, while for small values of 7, it manifests itself 
very clearly. 

The strength of the singularity is determined by the constant C, whose value 
depends on the input used. In particular, if the scattering length a is varied, 
while the absorptive part A(s) is kept fixed, the size of C changes. We may 
search for the value of a at which C vanishes. Although the singularity does not 
disappear entirely even then, it now only manifests itself in the derivatives of 
the function (for the solution to become analytic at s , we would need to also 
adapt the input for A(s)). In view of the fact that our solutions are inherently 
fuzzy, because the values of the input are subject to experimental uncertainties, 
we consider solutions with C ~ or 7 ~ 1 as physically acceptable and refer to 
these as solutions without cusp. 

The search for solutions without cusp can be implemented as follows. Instead 
of fixing a, constructing solutions in the class of functions with a cusp and then 
determining the value of a at which the cusp disappears, we may simply consider 
parametrizations that do not contain a cusp, treating the scattering length a as 
a free parameter, on the same footing as the set c±, . . . , c n used to parametrize 
the phase shift and minimizing the difference between the left and right hand 
sides of eq. Q6.1] ). We have verified that if a solution without cusp does exist, this 
procedure indeed finds it: Allowing for the presence of cusps does not lead to a 
better minimum. 

The net result of this discussion is that the scattering length a must match 
the input for A(s) - it does not represent an independent parameter. When 
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range of so 


range of 8q 


range of b\ 


m 


I 

II 
III 

IV 


1 < v^o < 115 
0.86 < < 1 
0.78 < v^< 0.86 
0.28 < v^< 0.78 


7T < 5° < | 7T 
\ 7T < 5° < 7T 
< 5° < \ 71 
< 5° < § 7T 


\ TV < 5{ < 7T 
| 7T < < 7T 
| 7T < $[ < 7T 
< 5\ < \ 7t 


2 
1 

-1 



Table 1: Multiplicity of solutions in the coupled channel case. The multiplicity 
index m is the number of free parameters occurring in the solutions of the Roy 
equations, if the matching point So is in the interval indicated (in GeV units). 
Also displayed is the variation of the physical phases 5$ and 5\ on that interval. 

solving the Roy equations, we can at the same time also determine the value 
of a that belongs to a given input for the high energy absorptive part. The 
conclusion remains valid even if the matching point is above the first inelastic 
threshold, provided the elasticity parameter rj is known and sufficiently smooth 
at the matching point. For a thorough analysis of the issue, we refer to [|31|j . 

6.3 Uniqueness in the mult i- channel case 

In the multichannel case, we need to determine three functions 5q, 5\ and 5$ for 
a given input {a®, a^, //, r/f , d{}. The multiplicity index m of the infinitesimal 
neighbourhood of a given solution is displayed in table [TJ PH , for various values 
of the matching point Sq. The table contains the following information. In the 
situations indicated with the labels I and II, the infinitesimal neighbourhood of 
a given solution contains a family of solutions, characterized by 2 and 1 free 
parameters, respectively. In case III, the solution is unique in the sense that the 
neighbourhood does not contain further solutions, while in case IV a solution only 
exists if the input is subject to a constraint (m = —1, compare paragraph 3 in 
section 6.1). In order to uniquely characterize the solution in case I, for instance, 
we thus need to fix two more parameters - in addition to the input - say the 
position of the p resonance and its width, or the position of the p resonance and 
the value of s where the 1 = phase passes through tt/2, and similarly for II. In 
the following, we stick to case III, where the solution is unique for a given input. 
As discussed above, each of the three partial waves will in general develop a cusp 
at the matching point s , unless some of the input parameters take special values. 

The situation encountered in practice is the following. Let 0.1 < a[| < 0.6, 
and let //, T]j and d\ be fixed as well. For an arbitrary value of the scattering 
length do, the solution in general develops a strong cusp in the P-wave. This 
cusp can be removed by tuning Oq — > a^, using for instance the method described 
in the single channel case above. Remarkably, it turns out that the solutions 
so obtained are nearly free of cusps in the two S-waves as well. The problem 
manifests itself almost exclusively in the P-wave, because our matching point is 



17 



rather close to the mass of the p, where the imaginary part shows a pronounced 
peak. If ag is chosen to slightly differ from the optimal value Oq, a cusp in the 
P-wave is clearly seen. We thus obtain a relation between the scattering lengths 
dp and Oq. This is how the so-called universal curve, discovered a long time ago 
[Pj] , shows up in our framework. We will discuss the properties of this curve in 
detail below. 

In principle, we might try to also fix a[] with this method, requiring that there 
be no cusp in one of the two S-waves. The cusps in these are very weak, however 
- the procedure does not allow us to accurately pin down the second scattering 
length. The choice a° = —0.2, for instance, still leads to a fully acceptable 
solution. On the other hand, we did not find a solution in the class of smooth 
functions for ct[] = —0.5. This shows that the analyticity properties that are not 
encoded in the Roy integral equations ( |5.1|) do constrain the range of admissible 
values for afj, but since that range is very large, the constraint is not of immediate 
interest, and we do not consider the matter further. In our numerical work, we 
consider values in the range 0.15 < a° < 0.30 and use the center of this interval, 
a[] = 0.225, as our reference point. 



7 Experimental input 

In this section, we describe the experimental input used for the elasticity below 
the matching point at y/so = 0.8 GeV and for the imaginary parts of the S- and 
P-waves in the energy interval between y/so and yfs~2 = 2 GeV. The references 
are listed in [fi5[|-||59| and for an overview, we refer to 0, |BP|| . The evaluation of 
the contributions from the higher partial waves and from the asymptotic region 
(s > s 2 ) is discussed in detail in appendix ||. 

7.1 Elasticity below the matching point 

The Roy equations allow us to determine the phase shifts of the S- and P- 
waves only if - on the interval between threshold and the matching point - 
the corresponding elasticity parameters rj^s), i]l(s) and tIq(s) are known. On 
kinematic grounds, the transition 2n — > An is the only inelastic channel open 
below our matching point, y/so = 0.8 GeV. The threshold for this reaction is 
at E = 4M n ~ 0.56 GeV, but phase space strongly suppresses the transition at 
low energies - a significant inelasticity only sets in above the matching point. In 
particular, the transition tttt — > KK, which occurs for E > 2Mk — 0.99 GeV, 
does generate a well-known, pronounced structure in the elasticity parameters of 
the waves with / = 0, 1. Below the matching point, however, we may neglect the 
inelastic reactions altogether and set 

r)o(s) = vl(s) = v 2 o(s) = 1 , < 0.8 GeV . 
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We add a remark concerning the effects generated by the inelastic reaction 
2tt — > 4tt, which are analyzed in ref. [[57). In one of the phase shift analyses given 
there (solution A), the inelasticity 1 — rj\(s) reaches values of order 4%, already 
in the region of the p -resonance. The effect is unphysical - it arises because the 
parametrization used does not account for the strong phase space suppression at 
the 47r threshold^. For the purpose of the analysis performed in ref. [[57]], which 



focuses on the region above 1 GeV, this is immaterial, but in our context, it 
matters: We have solved the Roy equations also with that representation for the 
elasticities. The result shows significant distortions, in particular in the P-wave. 

7.2 Input for the 7 = 0, 1 channels 

The experimental information on the irn phase shifts in the intermediate energy 
region comes mainly from the reaction ttN — > ttttN. A rather involved analysis is 
necessary to extract the 7nr phase shifts from the raw data, and several different 
representations for the phases and elasticities are available in the literature. The 
main source of experimental information is still the old measurement of the re- 
action ir~p — > 7T~7T + n by the CERN-Munich (CM) collaboration [fi~9[| , but there 
are also older, statistically less precise data, for instance from Saclay |45| and 



Berkeley |[48|| , as well as newer ones, such as the data of the CERN-Cracow- 



Munich collaboration concerning pion production on polarized protons [|54]] and 
those on the reaction n~p — > 7r°7r°n, obtained recently by the E852 collabora- 



tion at Brookhaven fSPfl . For a detailed discussion of the available experimental 
information, we refer to || |57|, |60[| . 

For our purposes, energy-dependent analyses are most convenient, because 
these yield analytic expressions for the imaginary parts, so that the relevant 
integrals can readily be worked out. To illustrate the differences between these 
analyses, we plot the corresponding imaginary parts in fig. E3, both for the 1 = 



S-wave and for the P-wave. The representations of refs. |47], |55|, |57| do not extend 
to 2 GeV, but they do cover the range between 0.8 and 1.7 GeV. Unitarity ensures 
that the contributions generated by the imaginary parts of the S- and P-waves 
in the region between 1.7 and 2 GeV are very small, so that we may use these 
representations also there without introducing a significant error. For the P-wave, 
the differences between the various parametrizations are not dramatic, but for 
the 1 = S-wave, they are quite substantial. Despite these differences, the result 
obtained for the dispersive integrals are similar, at least in the range where we 
are solving the Roy equations. This can be seen in fig. |], where we plot the value 
of the dispersion integral f$, defined in eq. (|5.2j ). The only visible difference is 



between parametrization B of ref. |57j and the others. In order of magnitude, the 
effect is comparable to the one occurring if the scattering length a{] is shifted by 
0.01. It arises from the difference in the behaviour of the .S-wave imaginary part 

4 We thank Wolfgang Ochs for this remark. 
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E(GeV) 



Figure 3: Comparison of the different input we used for the imaginary parts of 
the 1 = and 1 = 1 lowest partial waves above the matching point at 0.8 GeV. 



in the region between 1 and 1.5 GeV. The phase shift analysis of Protopopescu et 
al. JiS| does not cover that region, as it only extends to 1.15 GeV, but those of Au, 
Morgan and Pennington |55| as well as Bugg, Sarantsev and Zou |57j do. Both 
of these include, aside from the CM data, additional experimental information, 
not included in the analysis of Hyams et al. [47J. 

In the following, we rely on the representation of Au et al. pHI for the 5-wave 
and the one of Hyams et al. ^7| for the P-wave (the analysis of Au et al. does not 
include the P-wave). We have verified that, using [^7| also for the S'-wave would 
not change our results below the matching point, beyond the uncertainties to be 
attached to the solutions, anyway. On the other hand, Au et al. |55| yield a more 
consistent picture above the matching point - for this reason we stick to that 



analysis. More precisely, we use the solution denoted by ^(Etkin) in ref. [|55 
table I. That solution contains a narrow resonance in the 1 GeV region, which 
does not occur in the other phase shift analyses. In our opinion, the extra state 
is an artefact of the representation used: A close look reveals that the occurrence 
of this state hinges on small details of the K- matrix representation. In fact, the 
resonance disappears if two of the i^-matrix coefficients are slightly modified, for 
instance with (-c? 2 , -c° 22 ) = (3.1401,2.8447) -> (3.2019,2.6023). 
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Figure 4: Comparison of the results obtained for the dispersion integral /° with 
the various imaginary parts shown in fig. |3|. 



7.3 Phase of the P-wave from e 



7T ' 7T 



and r 



7T 7T° V T 



For the P-wave, the data on the processes e + e~ — > tt + tt~ and r — > 7r~7r° z/ T yield 
very useful, independent information. The corresponding transition amplitude is 
proportional to the pion form factor F e m (s) of the electromagnetic current and 
to the form factor F v (s) of the charged vector current, respectively. The data 
provide a measurement of the quantities |F e . m .(s)| and |-FV( S )| in the time-like 
region, s > AM 2 . 

In the isospin limit, the two form factors coincide: The currents only differ by 
an isoscalar operator that carries odd G-parity, so that the pion matrix elements 
thereof vanish. While the isospin breaking effects in |F v (s)| are very small, p — uo 
interference does produce a pronounced structure in the electromagnetic form 
factor. The tu-resonance generates a second sheet pole in the isoscalar matrix 
elements, at s = (M u — i^T^) 2 . The residue of the pole is small, of order 
0(rrid — m u ,e 2 ), but in view of the small width of the u, the denominator also 
nearly vanishes for s = M 2 . Moreover, the pole associated with the exchange of a 
p occurs in the immediate vicinity of this point, so that the transition amplitude 
involves a sum of two contributions that rapidly change with s, both in magnitude 
and phase. Since the interference phenomenon is well understood, it can be 
corrected for. When this is done, the data on the two processes e + e~ — > 7r + 7r~ 
and r — > 7t~7t v are in remarkably good agreement (for a review, see |61], 62|). 
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We denote the phase of the vector form factor by 4>(s), 

F v (s) = \F v (s)\e i ^ . 

In the elastic region 4M^ < s < 1GM%, the final state interaction exclusively 
involves nir scattering, so that the Watson theorem implies that the phase </>(s) 
coincides with the P-wave phase shift, 

(P(s) = 6{(s) , AMl<s<lQMl . 

In fact, phase space suppresses the inelastic channels also in this case - the 
available data on the decay channel r — > 4irv T show that, for E < 0.9 GeV, the 
inelasticity is below 1%, so that the phase of the form factor must agree with the 
P-wave phase shift, to high accuracy |63[ . 



In the region where the singularity generated by p -exchange dominates, in 
particular also in the vicinity of our matching point, the form factor is well rep- 
resented by a resonance term and a slowly varying background. Quite a few such 
representations may be found in the recent literature. Since the uncertainties 
in the data (statistical as well as systematic) are small, these parametrizations 
agree quite well. In the following, we use the Gounaris-Sakurai representation of 
ref. []64| as a reference point. That representation involves a linear superposition 



of three resonance terms, associated with p(770), p(1450) and p(1700). We have 
investigated the uncertainties to be attached to this representation by (a) com- 
paring the magnitude of the form factor with the available dataQ, (b) comparing 
it with other parametrizations, (c) varying the resonance parameters in the range 
quoted in ref. and (d) using the fact that analyticity imposes a strong cor- 



relation between the phase of the form factor and its magnitude. On the basis 
of this analysis, we conclude that the e + e~ and r data determine the phase of 
the P-wave at 0.8 GeV to within an uncertainty of ±2°. A detailed comparison 
between the phase of the form factor and the solution of the Roy equations for 



the P-wave will be given in section 12.2 



7.4 Phases at the matching point 

In the framework of our analysis, the input used for s > Sq enters in two ways: 
(i) it specifies the value of the three phases at the matching point and (ii) it 
determines the contributions to the Roy equation integrals from the region above 
that point. Qualitatively, we are dealing with a boundary value problem: At 
threshold, the phases vanish, while at the matching point, they are specified by 
the input. The solution of the Roy equations then yields the proper interpolation 
between these boundary values. The behaviour of the imaginary parts above the 
matching point is less important than the boundary values, because it only affects 
the slope and the curvature of the solution. 

5 We are indebted to Simon Eidclman and Fred Jegerlehncr for providing us with these. 
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s°0 

81.7 ± 3.9 
90.4 ± 3.6 
85.7 ± 2.9 



Si 



si 



105.2 ± 1.0 
115.2 ± 1.2 
116.0 ± 1.8 



23.4 ± 4.0 
24.8 ± 3.8 
30.3 ± 3.4 



reference 



50| s-channel moments 
50|| t-channel moments 



81.6 ± 4.0 108.1 ± 1.4 26.5 ± 4.2 



table VI 
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, 47 



f57J] solution A 



57] solution B 



55f| solution K] 



1 55[| solution Ki(Etkin) 



Table 2: Value of the phases 5$ and 5j at 0.8 GeV. The first three rows stem from 
analyses of the data at a fixed value of the energy ("energy independent"), while 
the remaining entries are obtained from a fit to the data that relies on an explicit 
parametrization of the energy dependence ("energy dependent analysis"). 



We now discuss the available information for the phases 5$ and S\ at the 
matching point. The values obtained from the high energy, high statistics ttN — > 
7nrN experiments are collected in table ||. In those cases where the published 
numbers do not directly apply at 0.8 GeV, we have used a quadratic interpolation 
between the three values of the energy closest to this one. The errors given in 
the third column are obtained by adding those from the first two columns in 
quadrature. For the energy dependent entries, the error analysis is more involved 
- only ref. |^8| explicitly quotes an error. The scatter seen in the table partly 



arises from the fact that different methods of analysis are used. The corresponding 
systematic uncertainties are not covered by the error bars quoted in the individual 
phase shift analyses: Taken at face value, the numbers listed in the table are 
contradictory, particularly in the case of the P-wave. For a thorough discussion 



of the experimental discrepancies, we refer to |60| . 

As discussed above, both the statistical and the systematic uncertainties of 
the e + e~ and r data are considerably smaller. They constrain the phase of the 
P-wave at 0.8 GeV to a narrow range, centered around the value Sl(so) = 108.9° 
obtained with the Gounaris-Sakurai representation of the form factor in ref. : 

5{(s ) = 108.9° ±2° . (7.1) 

The comparison with the numbers listed in the second column of the table shows 
that this value is within the range of the results obtained from ttN — > ttttN. 

Unfortunately, the e + e _ and r data only concern the P-wave. To pin down 
the 1 = S'-wave, we observe that the overall phase of the scattering amplitude 
drops out when considering the difference 5\ — 5®, so that one of the sources 
of systematic error is absent. Indeed, the third column in the table shows that 
the outcome of the various analyses is consistent with the assumption that the 
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fluctuations seen are of statistical origin. The statistical average of the energy 
independent analyses yields £i( s o) — <5q (s ) = 26.6° ± 3.7°, with x 2 = 2 for 2 
degrees of freedom (as the numbers are based on the same data, we have inflated 
the error bar - the number given is the mean error of the three data points). The 
remaining entries in the table neatly confirm this result. Combining it with the 
one in the fourth row, which is based on independent data, we finally arrive at 

<*i(*o) -$o(so) = 26.6° ±2.8° . (7.2) 

Since the value for 5{ comes from the data on the form factor, while the one for 
the difference 5\ — 5$ is based on the reaction ttN — > tcttN, these numbers are in- 
dependent, so that it is legitimate to combine them. Adding errors quadratically, 
we obtain 

5° {s ) = 82.3° ±3.4° . (7.3) 

In the following, we rely on the two values for the phases at the matching 
point given in eqs. (|7. 1| ) and (|7.3| ). We emphasize that the ttN — > ttitN data 
are consistent with these - in fact, the result of the energy-dependent analysis 
quoted in the fourth row of the table is in nearly perfect agreement with the 
above numbers. We are exploiting the fact that the e + e~ and r data strongly 
constrain the behaviour of the P-wave in the region of the p, thus reducing the 
uncertainties in the value of 5\ at the matching point. 

For the principal value integrals to exist, we need to continuously connect 
the values of the imaginary parts calculated from the phases at the matching 
point with those of the phase shift representation we wish to use. This can be 
done, either by slightly modifying the parameters occurring in the representation 
in question or with a suitable interpolation of the phases between the matching 
point and KK threshold. We have checked that our results do not depend on 
how that is done, as long as the interpolation is smooth. Note that, for the 
representation Ifi(Etkin) []55j - our reference input for the imaginary part of the 
1 = S-wave - an interpolation is not needed: The last row of table shows 
that, at the matching point, this representation nearly coincides with the central 
value in eq. ( |7.3| ). 



7.5 Input for the 1 = 2 channel 

The uncertainties in this channel are rather large. The current experimental 
situation is summarized in fig. [5], where we show the data points from the two main 



experiments |51], 53], and five different parametrizations that we will use as input. 



The central one is our best fit to the data of the Amsterdam-CERN-Munich 
collaboration (ACM) [B3J solution B (which we call from now on ACM(B)) with 
a parametrization a la Schenk [|65|] . To cover the rather wide scatter of the data, 
we have varied the input in this channel, using the five curves shown in the figure, 
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Figure 5: Different data sets for the S-wave in the 1 = 2 channel and curves 
that we have used as input in the Roy equation analysis. 

together with t]q = 1 (note that for the Roy equation analysis, only the value of 
the scattering length and the behaviour of the imaginary part above 0.8 GeV 
matter) . 

8 Numerical solutions 

In the preceding section, the input required to evaluate the r.h.s. of our system 
of equations was discussed in detail. In the present section, we describe the 
numerical method used to solve this system and illustrate the outcome with an 
example. 

8.1 Method used to find solutions 

We search for solutions of the Roy equations by numerically minimizing the square 
of the difference between the left and right hand sides of eq. (|5.1j) in the region 
between threshold and 0.8 GeV. As we are neglecting the inelasticity in this 
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region, the real and imaginary parts of t\(s) are determined by a single real 
function, the phase Sj(s). In principle, the minimization should be performed 
over the whole space of physically acceptable functions {Sq(s), Sl(s), 5q(s)}, but 
for obvious practical reasons we restrict ourselves to functions described by a 
simple parametrization. We will use the one proposed by Schenk some time ago 
|65|1 , allowing for an additional parameter in the polynomial part: 



tan 5j 



AM 2 r > f AM' 

^ 1V1 tt 21 ( tJ ■ nl 9. , nf J i^r fill ^ IVI l 



Q 



2 -si 



{Ai + Btf + C^ + Dtf}(^^\ , (8.1 



The first term represents the scattering length, while the second is related to the 
effective range: 

< = A \ > b l = B l + A l ~ ]^ bo ■ (8-2) 

In each channel, one of the five parameters is fixed in order to ensure the proper 
value of the phase at s . Moreover the S'-wave scattering lengths a° and a 2 , are 
identified with the two constants that specify the subtraction polynomials in the 
Roy equations. As discussed in sect. ^, we need to tune the value of in order 
to avoid cusps. Treating this parameter on the same footing as the others, we 
are dealing altogether with 15 — 3 — 1 = 11 free variables, to be determined by 
a minimization procedure. Our choice of so ensures that the solution is unique, 
and therefore the method is safe: The choice of a bad parametrization would 
manifest itself in a failure of the minimization method - the minimum would not 
yield a decent solution. 

The square of the difference between the left and right hand sides of the Roy 
equations is calculated at 22 points between threshold and s for each of the three 
waves, so that the sum of squares (A| oy ) contains 66 terms. The minimization 
of the function (A^ oy ) over 11 parameters can be handled by standard numerical 
routines Our procedure does generate decent solutions: The differences be- 



tween the left and right hand sides of the Roy equations are not visible on our 
plots - they are typically of order 10~ 3 . The equations could be solved even more 
accurately by allowing for more degrees of freedom in the parametrization of the 
phases, but, in view of the uncertainties in the input, the accuracy reached is per- 
fectly sufficient. Note also that the exact solution corresponding to a given input 
contains cusps. We have checked that these are too small to matter: Enlarging 
the space of functions on which the minimum is searched by explicitly allowing 
for such cusps in the parametrization of the phases, we find that the solutions 
remain practically the same. 



8.2 Illustration of the solutions 

To illustrate various features of our numerical solutions, we freeze for a moment 
all the inputs and analyze the properties of the specific solution we then get. 
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E(GeV) 

Figure 6: Numerical solution of the Roy equations for a° = 0.225, = —0.0371 
(the value of a[j corresponds to the center of the range considered while the one 
of <2q results if the input used for Imt^ is taken from the central curve in fig. |5|). 
The arrow indicates the limit of validity of the Roy equations. 



The input for the imaginary parts above So is the following: For the 1 = wave, 
we use the parametrization labelled K\ (Etkin) of Au et al. |5"5| . In the case of 
the 1 = 1 wave, we rely on the energy-dependent analysis of Hyams et al. [47], 
smoothly modified between s and 4M|- to match the value <5|(so) = 108.9°. For 
the 1 = 2 wave, we take the central curve in fig. [|. The driving terms are specified 
in eq. (|4.1|) . Moreover we fix a[j = 0.225. With this input, the minimization leads 
to Oq = —0.0371 and the Schenk parameters take the values listed in table [|, in 
units of M n . 

The plot in fig. ^| shows that the numerical solution is indeed very good: Below 
Sq, it is not possible to distinguish the two curves representing the right and left 
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Table 3: Schenk parameters of the solution shown in fig. ^. 
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hand sides of eq. flS.ip . For this solution we found as a minimum A^ oy = 2.1-10~ 5 , 
which corresponds to an average difference between the right and left hand sides 
of about 6 • 1CT 4 . 

Having solved the Roy equations in the low-energy region, we now have a 
representation for the imaginary parts of the three lowest partial waves from 
threshold up to &2- Since the driving terms account for all remaining contri- 
butions, we can then calculate the Roy representation for the real parts from 
threshold up to 1.15 GeV (full lines in fig.|]). On the same plot, above s , we also 
show the real part of the partial wave representation that we used as an input for 
the imaginary parts (dashed lines). The comparison shows that the input we are 
using is well compatible with the Roy equations (we should stress at this point 
that in none of the phase-shift analyses which we are using as input the Roy 
equations have been used). 



9 Universal band 

As we have discussed in the preceding sections, for a given value of a° and fixed 
input, the Roy equations admit a solution without cusp only for a single value 
of a 2 ,. By varying the input value of a®, the Roy equations define a function 



a 2 , = -F(a[]) that is known in the literature as the "universal curve" f4"4" |. The 
experimental uncertainties in the input above 0.8 GeV convert this curve into 
a band. The universal band is the area in the (oq, o,q) plane that is allowed 
by the constraints given by the 7T7r-scattering data above 0.8 GeV and the Roy 
equations. In this section we give a more precise definition of our universal band, 
and calculate it accordingly. 

We first point out that the universal curve a 2 , = F(a^) depends rather mildly 
on the input in the 1 = and 1 = 1 channel (a more quantitative statement 
concerning this dependence is given below). For this reason, we only consider 
the uncertainties in the input for the 1 = 2 channel. The available data in this 
channel are shown in fig. [|, together with five different curves that we have used as 
input. For each one of these, we obtain a universal curve, which nearly represents 
a straight line in the (a{j, a 2 ,) plane. The resulting five lines are shown in fig. [7[ 
The central one is well represented by the following second degree polynomial: 

a 2 Q = -0.0849 + 0.232 a° - 0.0865 (a° ) 2 . (9.1) 

The analogous representations for the top and bottom lines read: 



a 2 Q = -0.0774 + 0.240 a° - 0.0881 (a°) 2 , 

a 2 = -0.0922 + 0.225 a° -0.0847 (a°) 2 . (9.2) 

The region between these two solid lines is our universal band. It is difficult 
to make a precise statement in probabilistic terms of how unlikely it is that 
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Figure 7: Universal band. The five lines correspond to the five different curves 
shown in fig. [5| (the top line, for instance, results if the input for Im£g in the 
region above 0.8 GeV is taken from the top curve in that figure). So marks our 
reference point: = 0.225, Oq = —0.0371. The bar attached to it indicates the 
uncertainty in Oq due to the one in the phase Sq at the matching point - the most 
important remaining source of error if the input for Ihi£q is held fixed. 

the physical values of the two scattering lengths are outside this band. With 
our rather generous choice of the two extreme curves, we consider it fair to say 
that the experimental information above the matching point essentially excludes 
such values. In fact, we will argue below that the theoretical constraints arising 
from the consistency of the Roy equations above the matching point restrict the 
admissible region even further. 

We now turn to the dependence of the universal curve = on the 

input in the 1 = and 1 = 1 channels, keeping the one for / = 2 fixed. Changes 
in the input above 2M^ are practically invisible at threshold: If we keep the 
phase shifts at the matching point fixed, the three different available inputs for 
the 1 = and 1 = 1 channels yield values of Oq that differ by less than one 
permille. The phase shifts at so are the only relevant factor here. Moreover, for 
the value of Oq, Sq(s ) is much more important than <^(s ): Shifts of 5i( s o) by 
±2° change the value of roughly by a permille, but a change by ±3.4° in <5q(s ) 
induces a shift of Aa^ = ±8.4 ■ 10 -4 , which amounts to two percent. Even so, 
this is much smaller than the width of the band, as can be seen in fig. [F[ 
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We have also varied a/so within the bounds 0.78 and 0.86 GeV and found that 
the dependence of the relation a 2 , = F(cLq) on s is rather weak. To exemplify, we 
mention that for the solution with a{] = 0.25 at the center of the universal band, 
a shift from ^/s^ = 0.8 GeV to 0.85 GeV changes a 2 , by 10~ 3 . 

10 Consistency 

It takes a good balancing of the various terms occurring in the Roy equations for 
the partial waves not to violate the unitarity limit. In the case of the .S-wave with 
1 = 0, for instance, the contribution to Ret[j that arises from the subtraction term 
/cq(s) is very large already at 1 GeV: The solution shown in fig. ^| corresponds to 
a° = 0.225 and a 2 = -0.0371, so that k%(s) = 2.7 for s = 1 GeV 2 . As the energy 
grows, the term increases and reaches fcg( s i) = 3-6 at the upper end of the region 
where our equations are valid, si = 68 M 2 . Unless the contributions from the 
dispersion integrals nearly compensate the subtraction term, the unitarity limit, 
I Re | < (2<t) _1 — I is violated. The example in fig. || demonstrates that we do 
find solutions for which such a cancellation takes place, with values of af], a 2 , that 
are within the universal band. 

It is striking that, above the matching point, this solution very closely fol- 
lows the real part of the input. In a restricted sense, this is necessary for the 
solution to be acceptable physically: The solution is obtained by identifying the 
imaginary part above the matching point with the one obtained from a particular 
representation of the partial waves. The Roy equations then determine the real 
part of the amplitude in the region below ^/s± = 1.15 GeV. If the result were very 
different from the real part of the particular representation used, we would have 
to conclude that this representation cannot properly describe the physics. This 
amounts to a consistency condition: Above the matching point, the Roy solution 
should not strongly deviate from the real part of the input. The condition can 
be met only if the cancellation discussed above takes place, but it is stronger. 
The example in fig. |6] demonstrates that there are solutions that obey the consis- 
tency condition remarkably well, indicating that our apparatus is indeed working 
properly. 

We will discuss the consistency condition on a quantitative level below. Before 
entering this discussion, we briefly comment on a different aspect of our frame- 
work: the stability of the solutions. The behaviour below 0.8 GeV is not sensitive 
to the uncertainties in the input used for the imaginary parts above 1 GeV. We 
can modify that part of the input quite substantially, and without changing any- 
thing else (not even below sq) still get a decent solution from threshold up to 
the limit of validity of our equations. Naturally, if we do not modify the Schenk 
parameters that define the phase below s , the Roy equations are not strictly 
obeyed, but the deviation from the true solution is quite small. The reason is 
that, if s is small, the kernels K\l (s, s') strongly suppress the contributions from 
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Figure 8: Solutions of the Roy equations for a® = 0.225 and two extreme values 
for Oq. The left figure corresponds to the point S2 in fig.[?|, while the one on the 
right shows the solution for Si. The arrows indicate the limit of validity of the 
Roy equations. 



the region where s' is large. The term Kqq(s, s'), for instance, has the following 
expansion for s' s: 

s') = I {ll, 2 - 10 S (4M 2 ) - (4M 2 ) 2 } -L + . 

The interval above 1 GeV only generates very small contributions to the integrals 
on the r.h.s. of the Roy equations, if these are evaluated in the region below the 
matching point. 

We now take up the consistency condition and first observe that, once a solu- 
tion has a consistent behaviour above the matching point, reasonable changes in 
the input above 1 GeV lead to solutions that also obey the consistency condition: 
It looks as if the Roy equations were almost trivially satisfied, behaving like an 
identity for E > 1 GeV. Is this consistent behaviour automatic, or does it depend 
crucially on part of the input ? 

The answer to this question can be found in fig. where we show two solutions 
obtained with the same value of 0° as in fig. ^, but different inputs for hut 2 ,: The 
solution on the left is obtained by using the top curve in fig. ^| instead of the 
central one (a 2 , = —0.0279 instead of a 2 , = —0.0371). The solution on the right 
corresponds to the bottom curve in fig. [5], where a 2 , = —0.0460. The figure clearly 
shows that the consistent picture which we have at the center of the universal 
band is almost completely lost if we go to the upper border of this band: It is by 
no means trivial that we at all find solutions for which the output is consistent 
with the input. 

The fact that the peaks and valleys seen in the solutions mimic those in the 
input can be understood on the basis of analyticity alone: The curvature above 
the matching point arises from the behaviour of the imaginary parts there. The 
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relevant term is the one from the principal value integral, 



Ret(s) = - 4 



7T J AMI 



1 /"-a 



ds 



, Imt(s') 
s' — s 



+ r(s) 



The remainder, r(s) contains the contributions associated with the subtraction 
polynomial, the left hand cut, the higher partial waves, as well as the asymptotic 
region. On the interval s < s < si, it varies only slowly and is well approximated 
by a first order polynomial in s. 

The representations of the partial wave amplitudes that we are using as an 
input are specified in terms of simple functions. In the vicinity of the region 
where we are comparing their real parts with the Roy solutions, these are ana- 
lytic in s, except for the cut along the positive real axis. Hence they also admit 
an approximate representation of the above form - the contributions from distant 
singularities are well approximated by a first order polynomial. Disregarding the 
interpolation needed to match the representation with the prescribed value of 
the phase at sq, their imaginary parts coincide with the one of the corresponding 
Roy solution above the matching point. The small differences occurring in the 
interpolation region and below the matching point do not generate an important 
difference in the curvature. We conclude that the difference between the Roy so- 
lution and the real part of the input must be linear in s, to a good approximation. 
Moreover, within the accuracy to which our solutions obey the Roy equations, 
the two expressions agree at the matching point, by construction. Accordingly, 
the relation can be written in the form 



We have checked that this relation indeed holds to sufficient accuracy, for all three 
partial waves. This does not yet explain why the solution follows the real part of 
the input, but shows that it must do so up to a term linear in s that vanishes at 
the matching point. In particular, if the difference between input and output is 
small at the upper end of validity of our equations, then analyticity ensures that 
the same is true in the entire region between the matching point and that energy 
(in this interval, s varies by about a factor of two). 

In view of the uncertainties attached to our input, we cannot require the 
Roy equations to be strictly satisfied also above the matching point. The band 
spanned by the two green lines in fig. |9| shows the region in the (a[j, a*) plane, 
where the solution for Retg(s) differs from the real part of the input by less 
than 0.05 (expressed in terms of the parameter (3 in eq. ( |10.1|) , this amounts to 
\0q\ < 0.07GeV~ 2 ). Likewise, the band spanned by the two blue lines represents 
the region where |Reto(s) Roy - Reto( s )in P ut| < 0.05, so that |/3 2 | < 0.07 GeV~ 2 . 
The corresponding band for the P-wave is much broader - in this channel, the 
consistency condition is rather weak and is met everywhere inside the universal 
band. We conclude that, in the lower half of the universal band, all three waves 



Ret(s) 



Roy 



R ^( S ) input + (S-S )P • 



(10.1) 
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Figure 9: Regions inside which the consistency condition is met. The band 
between the two blue lines is for the condition in the 1 = 2 channel, whereas 
the one between the two green lines is for the 1 = channel. The two red lines 
delimit the band inside which the Olsson sum rule is satisfied. The shaded area 
gives the intersection of the three bands. 



show a consistent behaviour, while for the upper quarter of the band, this is not 
the case (the situation at the upper border is shown on the left in fig. §). 

It is not difficult to understand why the consistency condition is strongest 
for the 1 = S-wave. In this connection, the most important term in the Roy 
equations is the one from the subtraction polynomial - the solution can satisfy the 
consistency condition only if the term proportional to s is nearly cancelled by a 
linear growth of the remaining contributions. The term generates the contribution 
(Po,Pl,Po) = (6,1,-3) x (2al - 5ag)/(72A^) to the coefficients that describe 
the difference between output and input for the three lowest partial waves. The 
subtraction polynomial thus contributes twice as much to 0q as to (3%, so that the 
consistency band for the 1 = 2 wave must be about twice as broad as the one for 
the 1 = wave, while the one for the P-wave must roughly be six times broader. 
At the qualitative level, these features are indeed born out in the figure, but we 
stress that the term from the subtraction polynomial is not the only one that 
matters - those arising from the integrals also depend on the values of a(j and a^. 
The two green lines correspond to a variation in ag by about ±0.004. Increasing 
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Oq by 0.004, the value of the subtraction term &o(si) decreases by 0.10. The fact 
that the lines correspond to a change in Ret^si) of only ±0.05 implies that the 
contributions from the integrals reduce the shift by a factor of 2. Also, if only 
the subtraction term were relevant, the consistency bands would be determined 
by the combination 2a 2 , — 5a 2 , and thus have a slope of |. Actually, these bands 
are roughly parallel to the universal band, whose slope is positive, but smaller by 
about a factor of 2. 



11 Olsson sum rule 

In the Roy equations, the imaginary parts above the matching point and the two 
subtraction constants a[j, a 2 , appear as independent quantities. The consistency 
condition interrelates the two in such a manner that the contributions from the 
integrals over the imaginary parts nearly cancel the one from the subtraction 
term. In fact, a relation of this type can be derived on general grounds. 

The fixed-t dispersion relation ( |2.4| ) contains two subtractions. In principle, 
one subtraction suffices, for the following reason. The t-channel 1 = 1 amplitude 

TM(s, t) = \ {2 T°(s, t) + 3 T\s, t) - 5 T 2 (s, t)} 
o 

does not receive a Pomeron contribution and thus only grows in proportion to 
s a p(*) f or s _> oo. The dispersion relation (2.4), however, does contain terms that 
grow linearly with s. For the relation to be consistent with Regge asymptotics, the 
contribution from the subtraction term must cancel the one from the dispersion 
integral^. At t = 0, this condition reduces to the Olsson sum rule, which relates 
the subtraction constants to an integral over the imaginary parts ||67|| : 

2 Ml [°° 2ImT°(s,0) + 3ImT 1 (s,0) -5ImT 2 (s,0) 



It is well known that this sum rule converges only slowly - the contributions from 
the asymptotic region cannot be neglected. We split the integral into four pieces, 

2 a° - 5 a 2 Q = SP + D + F + O as . 

The first term represents the contributions from the imaginary parts of the S- 
and P- waves in the region below 2 GeV, which are readily worked out, using 
our Roy solutions on the interval from threshold to 0.8 GeV and the input phase 
shifts on the remainder. The result is not very sensitive to the input used and is 
well approximated by a linear dependence on the scattering lengths, 

SP = 0.483 ± 0.011 + 1.13 (a° - 0.225) - 1.01 (a 2 + 0.0371) . 



6 In the case of the t-channel amplitudes with I — and 1 = 2, the dxed-t dispersion relation 
(2.4) does ensure the proper asymptotic behaviour. 
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The remainder is closely related to the moments 2^ introduced in appendix B.l: 
here, we are concerned with the case n = — 1. The term D describes the con- 
tribution from the imaginary part of the .D-waves, in the interval from threshold 
to 2 GeV. The relevant experimental information is discussed in appendix B.3, 
where we also explain how we estimate the uncertainties. The numerical result 
reads D = 0.061 ± 0.004, including the small, negative contribution from the 
1 = 2 P-wave. The bulk stems from the tensor meson / 2 (1275): In the nar- 
row width approximation, this contribution amounts to 0.063. For the analogous 
contribution due to the F-wave, we obtain F = 0.017 ± 0.002 (in narrow width 
approximation, the term generated by the ps(1690) yields 0.013). Those from 
the asymptotic region are dominated by the leading Regge trajectory - as noted 
above, the Pomeron does not contribute. Evaluating the asymptotic contribu- 
tions with the formulae given in appendix B.4, we obtain O as = 0.102 ± 0.017. 
Collecting terms, this yields 

2 a ° - 5 a 2 = 0.663 ± 0.021 + 1.13 (ag - 0.225) - 1.01 {al + 0.0371) . (11.2) 

The result corresponds to a band in the (afj, al) plane: 

a 2 = -0.044 ±0.005 + 0.218 (ag - 0.225) . (11.3) 

The band is spanned by the two red lines shown in fig. |9]. One of these nearly 
coincides with the lower border of the universal band, while the other runs near 
the center. The Olsson sum rule thus imposes roughly the same relation between 
al and al as the consistency condition. Note that the asymptotic contributions 
are numerically quite important here: The term O as amounts to a shift in al 
of —0.026 ± 0.004. The fact that - in the region where our solutions are inter- 
nally consistent - the sum rule is indeed obeyed, represents a good check on our 
asymptotics. 

The Olsson sum rule ensures the proper asymptotic behaviour of the ampli- 
tude only for t = 0. In order for the terms that grow linearly with s to cancel 
also for t ^ 0, the imaginary part of the P-wave must obey an entire family of 
sum rules. The matter is discussed in detail in appendix C.I, where we demon- 
strate that one of these offers a further, rather sensitive test of our framework. 
The relationship between the Roy equations and those proposed by Chew and 
Mandelstam |68| is described in appendix C.2, where we also comment on the 



asymptotic behaviour of the dispersion integrals that occur on the r.h.s. of the 
Roy equations for the S- and P-waves. 



12 Comparison with experimental data 

In our framework, the only free parameter is ag. Comparing our Roy equation 
solutions to data, we can determine the range of ag consistent with these, as well 
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as a corresponding range for Oq. This experimental determination of the two 
S'-wave scattering lengths is the final scope of the present analysis and the main 
subject of the present section. Data on the 7T7T amplitude are available in a rather 
wide range of energies (we do not indicate the upper limit in energy when this 
exceeds 1.15 GeV, the limit of validity of our equations): 

• K eA data for the combination 5° - 8\ (2M W <E< 0.37 GeV); 

• ACM and Losty et al. data for <?g (0.35 GeV < E); 

• Data on the vector form factor - according to the discussion in section 
\l.3j these can safely be converted into values for 5{ in the region of the p 
(0.5 < E < 0.9 GeV); 

• CERN-Munich, and Berkeley data in the channels with 1 = and 1=1 
(0.5 GeV < E); 

In the Roy equations, a° and exclusively enter through the subtraction poly- 
nomials, specified in eq. (|1.2|) . Those relevant for the S- waves contain a con- 
stant contribution given by the scattering length and a term proportional to 
(s — 4:M%) x (2a{j — 5<2q). In the 1 = wave, that term is larger than a[] from 
E ~ 0.5 GeV on. For the 1 = 2 wave, the linear term starts dominating over 
Oq even earlier. Since t\(s) vanishes at threshold, the corresponding subtrac- 
tion polynomial exclusively involves the linear term. This implies that, except 
in the vicinity of threshold, the behaviour of the solutions is sensitive only to 
the combination 2a[J — 5oq of scattering lengths - roughly the combination that 
characterizes the universal band. Accordingly, only data that reach down close 
to threshold give a direct handle to separately determine a{] and Oq. In fact, only 
those coming from K e4: decays meet this condition. 

There is another threshold in energy that is obviously relevant for our ap- 
proach: the matching point sq. We will make a clear distinction between data 
points below so and those at higher energies. The comparison to data above so 
can hardly yield any information on the scattering lengths, because the behaviour 
of our solutions at those energies very strongly depends on the input used for the 
imaginary parts: The uncertainties in the experimental input completely cover 
the dependence of the solutions on the scattering lengths - we will discuss this in 
detail below. Instead, we analyze the requirement that the solution is consistent 
with the input for s > s , in the sense discussed in section [II]. This condition 
turns out to be practically independent of the input used for the imaginary parts 
above Sq and does therefore yield a meaningful constraint on 2afj — 5oq. 



12.1 Data on 5q — 5{ from K e 4, and on Sq below 0.8 GeV 

Let us first consider the K e ± data. The comparison between our solutions and the 
high-statistic data of the Geneva-Saclay collaboration [^] is shown in fig. [T^, for 
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Figure 10: Comparison of our Roy solutions for different values of the scattering 
lengths with the data of the Geneva-Saclay collaboration, Rosselet et al. [|S9 



The full, dash-dotted and dashed lines correspond to the points So, S2 and S3 in 
fig- 0- 



various values of the scattering lengths. The figure confirms the simple intuition 
that these data are mainly sensitive to a{j. In accordance with previous analyses 
7l| , we find that they roughly constrain a[J to the range between 0.18 and 0.3. 

As for the low-energy data in the 1 = 2 channel, we should stress that this 
wave is quite strongly constrained once Sq(so) is fixed. Because of the absence 
of any structures between threshold and 0.8 GeV, once we fix £ 2 (so), the only 
freedom is in the way the phase approaches zero at threshold, i.e. in the value 
of - which depends on a[j. Fig. [IT] shows that, at fixed Sq(s ), even a sizeable 
change in <2q is barely visible in the 1 = 2 phase. The only important factor 
here is the value of the phase at the matching point: The comparison with the 
experimental data basically tells us which value of Sq(sq) is preferred. 

A quantitative statement can be made in terms of x 2 > an d in principle we 
could calculate three different x 2_va hies, based on the three sets of data shown 
in fig. |5|. Two of these, however, represent two different analyses of the same 
set of ttN — > 7T7riV data. Their difference is a clear sign of the presence of 
sizeable systematic errors. We have estimated the latter using the difference, 
point by point, between the two analyses A and B of ref. p31, and added this in 
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Figure 11: Comparison of our Roy solutions with the data on 5q obtained by 
the ACM collaboration |k| and by Losty et al. |HJ. The full, dash-dotted and 
dashed lines correspond to the points So, S2 and S3 in fig. |7|. 



quadrature to the statistical errors. As reference we have used the ACM(A) set 
of data, but have checked that interchanging it with the one of Losty et al. does 
not give significantly different results. The corresponding % 2 , combined with the 
one obtained from the K ei data, has a minimum Xmin = 5.1 (with 8 d.o.f.) at 
tig = 0.242, <2q = —0.0357. The contour corresponding to 68% confidence level 
(x 2 = Xmin + 2-3) is shown in fig. |T2|: The range 0.18 < a§ < 0.3 is dictated by 
the K ei data, whereas the 1 = 2 data exclude the upper border of the band. 



12.2 The p resonance. 

The input used at the matching point implies that the P-wave phase shift must 
pass through 90° somewhere between threshold and 0.8 GeV - the Roy equations 
determine the place where this happens and how rapidly the phase must grow 
with the energy there. The solutions turn out to be very stiff: Varying the 
values of a[j and within the universal band, and also varying the input for the 
imaginary parts above 0.8 GeV within the experimental uncertainties, we obtain 
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Figure 12: Range selected by the data below 0.8 GeV. The dashed line represents 
the 68% C.L. contour obtained by combining the Geneva-Saclay data on K ei 
decay with those from ACM (A) on 8q. 

the narrow band of solutions shown in fig. [13|. 

In this figure, the energy range only extends to 0.82 GeV, for the following 
reason: Our solutions move along the Argand circle only below the matching 
point. At higher energies, the real part of the partial wave calculated from the 
Roy equations does not exactly match the imaginary part used as an input: 
unless we correct the latter, the elasticity r]\ differs from unity, already before 
the inelastic channels start making a significant contribution. If the consistency 
condition is met well, the departure from unity is small, but it can become as 
large as 5% if we go to the extreme of the consistency region shown in fig.^. This 
means that it does not make much sense to extract the value of the phase without 
adjusting the imaginary part. The proper way to do this is to extend the interval 
on which the Roy equations are solved, but we did not carry this out. 

In the region 0.7 GeV < 0.82 GeV, the result closely follows the data of the 
CERN-Munich collaboration. Below 0.7 GeV, however, the data are in conflict 
with the outcome of our analysis: The five lowest data points are outside the 
range allowed by the Roy equations, a problem noted already in ref. ||. In our 
opinion, we are using a generous estimate of the uncertainties to be attached to 
our input. Note, in particular, that at those energies, the driving terms barely 
contribute. We conclude that the discrepancy between our result and the CERN- 
Munich phase shift analysis occurring on the left wing of the p is likely to be 
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Figure 13: P-wave phase shift. The band shows the result of our analysis, 
obtained by varying the input within its uncertainties, while the data points 
indicate the phase shift measured in the process ttN — > ttttN by the CERN- 
Munich collaboration. The full line represents the phase of the vector form factor 
(Gounaris-Sakurai fit of ref. p4|). 



attributed to an underestimate of the experimental errors. As discussed below, 
the comparison with the e + e~ and r decay data corroborates this conclusion. 

Concerning the resonance parameters, we first give the ranges of mass and 
width that follow if, in the vicinity of the resonance, the phase shift is approxi- 
mated by a Breit-Wigner formula^] 

Ms) Ml + iY p M p -s T p M p 

e _ hj2 _-r< tit „ ' 



M 2 p -iT p M p - s ' to 1 w M 2 p -s ' 

In this approximation, the mass of the resonance is the real value of the energy 
where the phase passes through 90° and the width may be determined from the 
value of the slope d5\/ds at resonance. The solutions contained in the band shown 

7 The difference between ± iM p T p and (M p ± §r p ) 2 is beyond the accuracy of that 
approximation. The second is obtained from the first with the substitution M p — » M p — | T 2 p , 
MpTp — > MpT p , which increases the value of M p by about 4 MeV. 
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in the figure correspond to the range M p = 774 ± 3 MeV and T p = 145 ± 7 MeV, 
to be compared with the average values obtained by the Particle Data Group, 
M p = 770.0 ± 0.8 MeV, T p = 150.7 ± 1.1 MeV [KJ . 



The only process independent property of the resonance is the position of the 
corresponding pole - the above numbers specify this position only approximately. 
To determine it more accurately, we first observe that the Roy equations yield 
a representation of the partial wave t\(s) on the first sheet, in terms of the 
imaginary parts along the real axis. The first sheet contains both a right and 
a left hand cut. We need to analytically continue the function from the upper 
rim of the right hand cut into the lower half plane (second sheet). The difference 
between the values obtained in this manner and those found by evaluating the 
Roy representation in the lower half plane is given by the analytic continuation 
of the imaginary part, 

lmt](s) = —r^ sin 2 <^(s) . 
a{s) 

On the first sheet, t\(s) does not have singularities. Hence a pole can only 
arise from the continuation of the imaginary part. Indeed, the function sin 2 5j(s) 
contains the term exp 2 % 5j(s), which has a pole below the real axis. The position 
is readily worked out with the explicit, algebraic parametrization of the phase 
that we are using. The result illustrates an observation made long ago JTT], [T^, |75|| : 
The pole mass is lower than the energy at which the phase goes through 90°, by 
about 10 MeV: For the band shown in the figure, the pole position varies in the 
range 

M p = 762.5 ± 2 MeV , T p = 142 ± 7MeV . 

The e + e~ and r data neatly confirm the conclusion reached above: The phase 
of the form factor is in perfect agreement with the behaviour of the P-wave that 
follows from the Roy equations, but differs from the data of the CERN-Munich 
phase shift analysis, particularly below 0.7 GeV. In our opinion, the information 
obtained about the behaviour on the left wing of the resonance on the basis of the 
reactions e + e~ — > 7r + 7r~ and r — > 7r~7r°z/ is more reliable than the one obtained 
from ttN — > tttcN. The fact that the Roy equations are in good agreement with 
the e + e~ and r data is very encouraging. 

In view of the clean determination of the P-wave phase shift through e + e~ 
and r experiments, we find it instructive to draw fixed x 2 -contours in the (ag, clq) 
plane. To do so, we first need to attach an error bar to the curve representing 
the phase shift. In section [7.4|, we estimated the uncertainty in 5j(so) at ±2° or 



±2%. As we go down in energy, the relative precision of the determination of 
the phase decreases: A generous estimate of the uncertainty at ^/s = 0.5 GeV is 
10% or ±0.6°. A smooth interpolation between these two values is our estimate 
of the experimental error bar (below that energy, the e + e~ and r data become 
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Figure 14: 68% C.L. contour obtained by combining all relevant low energy 
data: K e4 decay, ACM(A) data on 5$ below 0.8 GeV and results for 6{ extracted 
from the e + e~ and r data on the pion form factor. 

scarce and have sizeable uncertainties). To construct the x 2 we have compared 
our solutions to the experimentally determined phase shift at five points between 
0.5 and 0.75 GeV. Combining this x 2 with those from the data on K e4 decays 
and on 5q below 0.8 GeV, we obtain the 68% C.L. area drawn in fig. |14[ The 
minimum of the x 2 is now 5.4 (with 13 d.o.f.). The position of the minimum is 
barely shifted: It now occurs at a® = 0.240, = —0.0356. In other words, at the 
place where the x 2 °f the K e4 data on Sq — 8{ and those on 5q had a minimum, 
the x 2 relative to the data on the form factor is practically zero and also has a 
minimum. In view of the fact that the uncertainties in 8{ are very small, this 
is quite remarkable. The data on the P-wave do not change the position of 
the minimum, but shrink the ellipse along the width of the universal band. As 
expected, they do not reduce the range of allowed values of a[}. 



12.3 Data on the I = 5-wave below 0.8 GeV 

In fig. [15] we compare the 5-wave obtained from our Roy equation solutions with 
the available data: CERN-Munich j|7| and Berkeley |48|]. The band shown is a 
representation of the uncertainties in the solution, which have two main sources: 
the uncertainty in 5g(s ) and the one in 8q(s ) (width of the universal band). 
The central curve shows our reference solution a[j = 0.225, <2q = —0.0371. The 
uncertainties indicated do not account for the changes occurring if the value 
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E(GeV) 

Figure 15: Comparison between the Roy solution for the S"-wave and the phase 
shift analyses of the CERN-Munich (circles) and Berkeley (squares) collabora- 
tions. The band shows the uncertainties in the Roy solution, which are dominated 
by those in ctQ and Oq. 



a[] = 0.225 is modified. Changing this value within reasonable bounds, however, 
brings the solution out of the band only below 0.4 GeV, already far below the 
first data point. The figure shows good agreement with the data, especially so 
for the Berkeley data set. The CERN-Munich data set shows a certain structure, 
which does not occur in our solutions - in view of the uncertainties in the data, 
this difference does not represent a problem. 

Despite the positive picture which emerges from the comparison, we refrain 
from using these data to draw confidence-level contours in the (a[j, Oq) plane. 
The 5-wave phase shifts have been extracted simultaneously with the P-wave. 
As discussed in the preceding section, these are affected by systematic errors 
which are at least as large as the statistical ones. The same must be true for 
the data in the 1 = channel, so that a quantitative comparison with the Roy 
solutions is barely significant. 



12.4 Data above 0.8 GeV 



The Roy equations are valid up to y/si = 1.15 GeV. In fig. [16], we show three 
different solutions for the 1 = and 1 = 1 partial waves, in the region above 
the matching point. They are obtained by using three different inputs for the 
imaginary parts (note that the curves represent our solutions, not the real parts 
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Figure 16: Behaviour of the solutions above the matching point. The curves 
show the solutions obtained with three different inputs for the imaginary parts. 
The data points are taken from the energy independent analysis of the CERN- 



Munich data [47]. The 1 = S'-wave is shown in black, the 1=1 P-wave in 
blue. 



of the input). The figure shows that the differences are substantial, especially in 
the S-wave, despite the fact that, below a/so = 0.8 GeV, the three solutions are 
practically identical, for all three waves. Evidently, above the matching point, 
the Roy solutions are very sensitive to the input used for the imaginary parts. 

It is not difficult to understand why that is so. As discussed in detail in section 
10, the solutions follow the real parts of the representation that is used as input 
(see fig. H for the case of Au et al. - in the other two cases, the picture is similar). 
The real parts of the three representations differ considerably. Moreover, all of 
these are systematically lower than the "data points" in fig. O, which show the 



result of an energy independent analysis of the CERN-Munich data ||47]|. In view 



of this, it is not surprising that the three Roy solutions are quite different and 
that they are also systematically lower than the data points. 

We conclude that a comparison of the Roy solutions with the data in the region 
above the matching point does not yield reliable information about the values of 
the two S'-wave scattering lengths and we do therefore not show confidence-level 



44 



contours relative to data above 0.8 GeV. 



13 Allowed range for a§ and Oq 

The above discussion has made clear that we can rely only on two rather solid 
sources of experimental information to determine the two S"-wave scattering 
lengths: the data on K ei and those on the P-wave in the p region. The for- 
mer determine a range of allowed values for a[] while the latter yield a range for 
the combination 2a{j — 5<2q. The consistency condition and the Olsson sum rule 



impose further constraints. Figure 17 summarizes our findings: We have super- 



imposed the ellipse of fig. [H] with the lines that delimit the consistency bands for 
the two S'-waves, as well as those relevant for the Olsson sum rule. The allowed 
range for a° and is the intersection of the ellipse with the band where the Ols- 
son sum rule is obeyed within the estimated errors. In that region, the solutions 
also satisfy the consistency condition. 
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Figure 17: Intersection of the ellipse in fig. [H] (68% CL. relative to the data 
on K e 4 decay, on 5$ and on the form factor) with the bands allowed by the 
consistency condition in all the three channels and by the Olsson sum rule. 



We find it quite remarkable that the data on the shape of the p resonance, the 
consistency condition and the Olsson sum rule all show a preference for the lower 
part of the universal band. This gives us confidence that our conclusion on which 
region in the (a®, aft) plane is allowed by the present experimental information is 
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rather solid. Once the new data on K e ^ decays will become available, the allowed 
range in a[j will become much narrower, and we will have a very small ellipse. The 
prospects of making a real precision test of the predictions for the two S-wave 
scattering lengths in the near future, appear to be very good, in particular also 



in view of the pionium experiment under way at CERN [29 



The 7iN — > 7mN data do provide essential information concerning the input 
of our calculations, but, as discussed in sections |12.3| and |12.4j , they do not impose 
a firm constraint on the scattering lengths (incidentally, these data also prefer 
the lower half of the universal band). This is unfortunate, because the power of 
the Roy equations (unitarity, crossing symmetry and analyticity) is that of con- 
necting regions of very different energy scales. The behaviour of the two S-waves 
in the immediate vicinity of threshold is determined by the scattering lengths. In 
the combination 2a|] — 5oq, these also determine the linear growth of the subtrac- 
tion polynomial: As we discussed in detail in section 10, the large contribution 
from the polynomial must be compensated to a high degree of accuracy by the 
dispersive integrals. We therefore expect that a reanalysis of the irN — > 7nrN 
data based on the Roy equations would lead to a rather stringent constraint on 
the allowed region, as it would make full use of the information contained in these 
data - in our opinion, the existing phase shift analyses are a comparatively poor 
substitute. 



14 Threshold parameters 
14.1 S- and P- waves 



As shown in ref. |T4" |, the effective ranges of the S- and P- waves and the P- 
wave scattering length can be expressed in the form of sum rules, involving inte- 
grals over the imaginary parts of the scattering amplitude and the combination 
2oq — 5oq of S'-wave scattering lengths. The sum rules may be derived from the 
Roy representation by expanding the r.h.s. of eq. (|5.1|) in q 2 and reading off the 
coefficients according to eq. (|2.3|) . In the case of the S'-wave effective ranges, the 
expansion can be interchanged with the integration over the imaginary parts only 
after removing the threshold singularity. This can be done by supplementing the 
integrand with a term proportional to the derivative 



d 1 h(s) 



h{s) = {s-2Ml)^s{s-AMl) . 



ds V / S ( S -4M2) {s(s-AMl)Y 
In this notation, the sum rules may be written in the form: 

- 9M 2 (s - AMI) MOO + 5M 2 (s - AM 2 W ) Imt%(s) - f (a°) 2 h(s)} 
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'a° Th(s)+b° 0d , (14.1) 



nJ S2 {s(s-AMl)Y 

b ^-ek « - ^ + 1 jfc W^mW ^ ~ AMl) lmt ° ois) 

+ 9M*(s - AMI) Imt}(s) + M 2 (7s - AM 2 ) Imtg(s) - 3(a 2 ) 2 h(s)} 
ds 



nJ S2 {s{s-AMl)Y 



<YKs)+bl 



d > 



+ 9(3s - 4M^) lmt\(s) + 5(s - 4M^) Imtjj(s)} + a} d , 

6 ! ^ il C {.(. -*m?)}> {" 2(s " 4A& ImtS(s) + 9 ( 3 s3 " 12 

+ 48sM*-64M^) Imt}(s) + 5(s-4M2) 3 Imt2( s )} + 6} d . 

The integrals only involve the imaginary parts of the S- and P-waves and are 
cut off at s — S2- The contributions from higher energies, as well as those from 
the imaginary parts of the partial waves with i = 2, 3, . . . are contained in the 
constants &Q d , bl d , a\ d , b\ d . By construction, these represent derivatives of the 
driving terms at threshold, 

d° (s) = q %° d + O( q *), d\(s) = q 2 a\ d + q %\ d + 0(q e ), d*{s) = q%l d + 0{<f). 

The numerical values obtained within our framework are given in the upper 
half of table H, where we also show the numbers quoted in the compilation of 
Nagels et al. JnJ, which are based on the analysis of Basdevant, Froggatt and 
Petersen ||. In accordance with the literature, we use pion mass units. Since 
the relevant physical scale is of the order of 1 GeV, the numerical values rapidly 
decrease with the dimension of the quantity listed. The columns A - E indicate 
the following contributions to the totalQ 

A. Contribution from the subtraction term oc 2oq — 5<2q. 

B. Imaginary parts of the S- and P-waves on the interval AM 2 < s < sq. This 

contribution is evaluated with the Roy solutions described in the text. 

C. Imaginary parts of the S- and P-waves in the range sq < s < S2- Here, we 

are relying on the experimental information, discussed in section [7]. 

D. Imaginary parts of the higher partial waves in the range AM 2 < s < s%. 

These are calculated in the same manner as for the driving terms of the S- 
and P-waves (see appendix |B.3| ) . 



8 The numbers given for the total include the tiny additional contributions to 6q and that 
arise from the integrals over /i(s)(oq) 2 and h(s)(a,Q) 2 in the interval S2 < s < oo. Numerically, 
these amount to 6b% = -6.3-1Q- 4 M~ 2 and 5b 2 = -1.7-10" 5 M~ 2 . 
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A 


B 


C 


D 


E 


total 


Ai 


A 2 


ref. 


75 




units 




2.12 


.45 


-.03 


.02 


.00 


2.56 


±.02 


+ .28 
-.12 


2.5 ±0.3 


10- X M- 2 


bl 


-1.06 


.26 


.02 


.01 


.00 


-.77 


±.003 


+ .03 
-.07 


-.82 ± . 


08 




a\ 


3.53 


-.03 


.13 


-.01 


.01 


3.63 


±.02 


+ .29 
-.11 


3.8 ±0.2 


10- 2 M- 2 


b\ 




4.05 


1.39 


-.07 


.08 


5.45 


±.13 


+ .35 
-.44 








10- 3 M" 4 


al 




1.29 


.28 


.07 


.03 


1.67 


±.01 


+ .15 
-.06 


1.7 ± . 


3 


lO^M" 4 


bl 
a 2 2 




-3.48 
1.67 


-.04 
-.51 


.25 
.35 


.02 
.02 


-3.25 
1.53 


±.07 
±.07 


+ .34 

-.87 

+1.1 
-.45 


1.3 ±3 


10 -4 M -6 
10 -4 M -4 


b 2 




-3.10 


-.09 


.06 


.02 


-3.11 


±.07 


+ .41 
-.95 








10- 4 M- 6 


a\ 
b 1 

°3 




5.11 
-3.96 


.26 
-.01 


.05 
.01 


.01 
.01 


5.43 
-3.95 


±.1 
±.08 


+1.6 
-.72 
+.89 
-1.9 


6±2 


10" 5 M- 6 
10~ 5 M- 8 



Table 4: Threshold parameters of the S-, P-, D- and F-waves. The significance 
of the entries in columns A-E is specified in the text. The column Ai indicates 
the uncertainty due to the error bars in the experimental input at and above 0.8 
GeV, whereas A 2 shows the shifts occurring if a® and a^ are varied within the 
ellipse of fig. IH, according to eqs. (|14.2|) and (|14.4j). 



E. Asymptotic contributions, s > s 2 . These are evaluated with the representa- 



tion given in appendix B.4 



For the reasons discussed earlier, we use y/so = 0.8 GeV, ^/s2 = 2 GeV. The 
values quoted in columns A and B are obtained with our reference solution, 
dp = 0.225, Oq = —0.0371, which corresponds to the point So in fig. [7|. 

The table shows that the result for 6°, 6q, a\, b\ is dominated by the contri- 
butions from the subtraction term and from the imaginary parts of the S- and 
P-waves. The higher partial waves and the asymptotic region only yield tiny 
corrections. The sum D±E represents the contribution from the driving terms. 
In the evaluation of these terms, which is discussed in detail in appendix [B.5| , we 
have constrained the polynomial fit with the relevant derivatives at threshold, so 
that the numerical values of the four constants b^ d , b^ d , a\ d , b\ d are correctly 
reproduced by the corresponding terms in the representation ( fO ). 



The uncertainty given in column Ai of table [| only accounts for the noise 
seen in our evaluation for the specific values a[] = 0.225, Oq = —0.0371 (errors in 
columns B-E added up quadratically). The sensitivity to these two parameters 
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is well represented by linear relations of the formf] 



b° = 2.56 x lO^M- 2 { 1 + 3.2 Aa° -12.7 Aa 2 } , 
W 2 = -0.77 x 10-^M- 2 { 1 + 2.5 Aa° - 7.6 Aa 2 } , 
a}= 3.63 x 10- 2 M- 2 {1 + 2.3 Aa° -7.8 Aag } , 
6}= 5.45xlO" 3 M- 4 {1 + O.lAag -5.7Aa 2 }, 

with Aafj = a[] — 0.225, Aa 2 , = a 2 , + 0.0371. Using this representation, the la 
ellipse of fig. [H] can be translated into la ranges for the various quantities listed 
in the table - these are shown in column A 2 (since our reference point is not at 
the center of the ellipse, the ranges are asymmetric). 

The table neatly demonstrates that the two S-wave scattering lengths are 
the essential low energy parameters - the uncertainty in the result is due almost 



exclusively to the one in a°, a 2 ,- This is to be expected on general grounds [[76 



The integrals occurring in the above sum rules are rapidly convergent, so that 
only the behaviour of the partial waves in the threshold region matters. The 
uncertainties in the input used for the imaginary parts above the matching point 
only enter indirectly, through their effect on the S- and P-waves in the threshold 
region. We did not expect, however, that the effect would be as small as indicated 
in the table and add a few comments concerning this remarkable finding. 

In order to document the statement that the uncertainties which we are at- 
taching to the phenomenological input of our calculation (behaviour of the imagi- 
nary parts above the matching point, elasticity, driving terms) only have a minute 
effect on the result for the threshold parameters, we find it best to give the nu- 
merical size of this effect (column Ai of the table). We repeat that the numbers 
quoted there merely indicate the noise seen in our evaluation - we do not claim to 
describe the scattering amplitude to that accuracy. Isospin breaking, for instance, 
cannot be neglected at that level of precision. 

The reason why the threshold parameters are insensitive to the uncertainties 
of our input is the following. As discussed in detail in sections |6|-|9], the solutions 
of the Roy equations in general exhibit a cusp at the matching point. If the 
imaginary parts above 0.8 GeV and the value of afj are specified, there is a solution 
with physically acceptable behaviour in the vicinity of the matching point only if 
the parameter a 2 , is chosen properly. In other words, there is a strong correlation 
between the behaviour of the imaginary parts and the parameters a[], a 2 ,. As we 
are selecting a specific value for these parameters, we are in effect subjecting the 
imaginary parts to a constraint. For this reason, the uncertainties in the input 
can barely be seen in the output for the threshold parameters - the main effect 
is hidden in a®, a^. The correlation just described originates in the fact that one 



9 For 0.15 < <2q < 0.30 the representation holds inside the universal band to better than 4%. 
Similar relations also follow directly from the representation of the S- and P-waves given in 
appendix but since the threshold region does not carry particular weight when solving the 
Roy equations, these do not have the same accuracy. 
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of the two subtraction constants is superfluous: The combination 2 a° — 5 ag may 
be represented as a convergent dispersion integral over the imaginary part of the 
amplitude. 

The correlation is illustrated by the lines in fig. |7], which correspond to the 
specific parametrization of the input used for the imaginary part of the 1 = 2 
S'-wave shown in fig. [5[ As there is very little experimental information about 
the energy dependence of this partial wave, we have worked out the change in 
the Roy solutions that occurs if this energy dependence is modified above the 
matching point. The result for the threshold parameters turns out to be practi- 
cally unaffected. Also, we have varied the driving terms within the uncertainties 
given in section (|. Again, the response in the threshold parameters can barely 
be seen. 



14.2 D- and F- waves 

Similar sum rules also hold for the threshold parameters of the higher partial 
waves. The contributions from the imaginary parts of the S- and P-waves are 
obtained by expanding the kernels occurring in the Roy equations for the D- and 
F-waves around threshold. We write the result in the form 

+ 5(s-4M„ 2 )ImijS( s )} + < i , 
+ 5(>-4JM|)ImtS(»)} + ^ , 

457T J4M| S d [S — AM*) L 

+ (s-4M*)Imtl(s)} + al d , (14.3) 
bl = ~ f\ 4 , {2(s - AMI) lmt° Q (s) + 3(s - 12M 2 ) lmt\(s) 

15 7T J AM% S [S — AM£) L 

+ (s-AM 2 n )lmtl(s)} + bl d , 
-5(s-AM 2 w )lmt 2 (s)} + al d , 
-5(s-AM 2 )lmt 2 (s)} + bl d , 
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where a^, b^ d , . . . contain the contributions from s > S2 as well as those from 
the higher partial waves. The evaluation of these contributions, however, meets 
with problems that we need to discuss in some detail. 

First, we note that the definition of the driving terms in eq. ( |3.2|) is suitable 
only for the analysis of the S- and P-waves. For £ > 2, the functions d\(s) contain 
a branch cut at threshold, so that these quantities are complex. In order to solve 
the Roy equations for the P-waves, for instance, the contributions generated by 
their imaginary parts need to be isolated, using a different decomposition of the 
right hand side of these equations. As far as the scattering lengths and effective 
ranges are concerned, however, only the values of the functions d\(s) and their 
first derivatives at threshold are needed, which are real. 

A more subtle problem arises from the fact that the explicit form of the 
kernels occurring in the Roy equations for the higher partial waves depends on 
the choice of the partial wave projection. As discussed in detail in ref. fTFfl , the 



definition (|A.4j ) - which we used in our analysis of the S- and P-waves - does not 
automatically ensure that the threshold behaviour of the partial waves with I > 3 
starts with the power q 21 . The problem arises from the fact that the solution of 
the Roy equations leads to a crossing symmetric scattering amplitude only if the 
imaginary parts of the higher partial waves satisfy sum rules such as the one in 
eq. (|B.8| ). In particular, the expansion of the P-wave in powers of q in general 
starts with 

Ret\{s) = x\q A + a\q Q + b\q* + . . . 

For the fictitious term x\ to be absent, the imaginary parts of the higher partial 
waves must obey a sum rule. In fact, we have written down the relevant sum 
rule already: equation (TB.8|) . The derivation given in section |B.2| shows that this 



constraint ensures crossing symmetry of the terms occurring in the expansion of 
the scattering amplitude around threshold, up to and including contributions of 
0(g 4 ). The threshold expansion of the partial waves with i > 3 thus only starts 
at 0(q 6 ) if this condition holds - in particular x\ then vanishes. The sum rule 
that allows us to pin down the asymptotic contributions to the driving terms 
for the S- and P-waves thus at the same time also ensures the proper threshold 
behaviour of the P-waves. The absence of a term of 0(q e ) in the G- waves leads 
to a new constraint, which could be derived in the same manner, etc. Note that 
the contributions from the imaginary parts of the S- and P-waves are manifestly 
crossing symmetric - the constraints imposed by crossing symmetry exclusively 
concern the higher waves^. 

The P-wave scattering length occurs in the expansion of the amplitude around 
threshold among the contributions of 0(g 6 ), two powers of q beyond the term just 

10 The family of sum rules discussed in appendix |C.l| does not follow from crossing symmetry, 
but from an asymptotic condition that goes beyond the Roy equations. As shown there, those 
sum rules do tie the imaginary part of the P-wave to the higher partial waves. 
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discussed. In the numerical analysis, we thus need to make sure that the sum rule 
holds to high precision if we are to get a reliable value in this manner. For the 
effective range, the situation is even worse. This indicates that for the numerical 
analysis of the higher partial waves, the extension of the range of validity of the 
Roy equations achieved if the standard partial wave projection ( |A.2j ) is replaced 
by (|A.3|) generates considerable complications. 

For the evaluation of the threshold parameters, this extension is not needed 
- we may use the partial wave projection (|A.2| ), for which the problem discussed 
above does not occur. In particular, x\ then automatically vanishes, so that the 
evaluation of the scattering lengths and effective ranges does not pose special 
numerical problems. To evaluate those from the asymptotic region, we expand 
the fixed-t dispersion relation ( |2.4| ) in powers of t. The results obtained for 



aP = 0.225, a 2 G = -0.0371 are listed in the lower half of table |. 

The dependence on the S'-wave scattering lengths may again be represented 
(to better than 6% inside the universal band for 0.15 < < 0.30) with a set of 
linear relations: 



al = 


1.67 x 10 


bl = 


-3.25 x 10 


a 2 = 


1.53 x 10 


b 2 - 


-3.11 x 10 


a l = 


5.43 x 10 


b 1 - 


-3.95 x 10 



-4 



3 M- 4 {1 + 2.6 Aag-8.6 Aa 2 } , 
M~ e {l + 6.6Aa° -17 A^ } , 
4 M" 4 {1 + 14Aa°-25Aa2}, 
4 M" 6 {1 + 6.2 Aa° -11 Aa 2 } , 1 > 

5 M~ 6 {l + 5.5Aa° -8Aa2 } , 
5 M" 8 {1+ 8Aa° -8Aa 2 } . 



The sensitivity is more pronounced here than in the case of the threshold pa- 
rameters for the S- and P-waves. In particular, the linear representation for the 
.D-wave scattering length a 2 , only holds to a good approximation if and do 
not deviate too much from the central values. 



15 Values of the phase shifts at s = M\ 

A class of important physical processes where the tot phase shifts play a relevant 
role is that of kaon decays. Let us recall, for instance, that the phase of e' is 
given by the value of 5 2 — 5q + \ tt at s = M\. In this section, we give numerical 
values for the three phase shifts at the kaon mass as they come out from our Roy 
equation analysis, and show the explicit dependence on the two S'-wave scattering 
lengths. In this manner, an improved determination of the latter will immediately 
translate into a better knowledge of the phases at s = M\_. 

The decays K° — > nn and K + — > -nit concern slightly different values of the 
energy. In view of the fact that the CP-violating parameter e' manifests itself in 
the decays of the neutral kaons, we evaluate the phases at s — M^ . Note that, 
in addition to this difference in the masses, there are also isospin breaking effects 
in the relevant transition matrix elements. As far as the im phases are concerned, 
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Value at 

s = M| 


Ax 


A 2 




37.3 


±1.4 


+4.3 
-1.6 


s\ 


5.5 


±0.1 


+ .3 
-.13 


si 


-7.8 


±0.04 


+ .7 

-.8 


S° ~ SI 


45.2 


±1.3 


+4.5 
-1.6 



Table 5: Values of the phase shifts at s = M^o in degrees. The central value is 
obtained with our reference solution of the Roy equations, where a[j = 0.225, = 
—0.0371. The column A x indicates the uncertainty due to the error bars in the 
experimental input at and above 0.8 GeV, whereas A 2 shows the shifts occurring 
if a° and Oq are varied within the ellipse of fig. [14], according to eq. (|15.1| ). 



however, the isospin breaking effects due to — m u are tiny, because G-parity 
implies that these only occur at order (rrid — m u ) 2 . 

As in the preceding section, we give values at the reference point a{j = 0.225 
and Oq = —0.0371, and break down the errors into those due to the noise in 
our calculations and those due to the poorly known values of the two scattering 
lengths. The results are shown in table [|. Like for the threshold parameters, the 
two S'-wave scattering lengths are the main source of uncertainty. In the present 
case, the errors due to the uncertainties in our experimental input at 0.8 GeV 
are not negligible, but they amount to at most 4%. 

The dependence of the central values on the two scattering lengths is well 
described by the following polynomials: 

^(M| ) = 37.3° {l + 3.0Aa° - 8.5Aa^} , 
8\{M 2 K o) = 5.5° {l + 1.7Aa° - Q.7Aa 2 } , (15.1) 
5 2 {M 2 K o) = -7.8° {l + 1.9Aag - 13Aa'} , 
5° {M 2 K0 ) - 5 2 {M 2 KO ) = 45.2° {l + 2.8Aa° - 9.4Aa'} . 

Our results are in agreement with refs. |)0], [78|, [79| , but are more accurate. In 
the foreseeable future, the two .S'-wave scattering lengths will be pinned down to 
good precision, so that the above formulae will fix the phases to within remarkably 
small uncertainties. 



16 Comparison with earlier work 

The Roy equations were used to obtain information on the tctt scattering ampli- 
tudes, already in the early seventies. Most of the work done since then either 
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follows the method of Pennington and Protopopescu |3], |4]] or the one of Basde- 
vant, Froggatt and Petersen || ||. In the present section, we briefly compare 
these two approaches with ours. A review of the results obtained by means of 
the Roy equations is given in ref. [f§. 

To our knowledge, Pennington and Protopopescu H were the first to analyze 
7T7T scattering data using Roy's equations. In principle, the approach of these 
authors is similar to ours. In our language, they fixed the matching point at 
i/io = 0.48 GeV. As input data, they relied on the 7T7T production experiment of 
the Berkeley group |^8| , using the data of Baton et al. |^5| for the 1 = 2 channel 
(at the time they performed the analysis, the high-energy, high-statistics CERN- 
Munich data [f|7| were not yet available). The Roy equations then allowed them 
to extrapolate the S- and P-wave phases of Protopopescu et al. [48] to the region 
below 0.48 GeV. Comparing the Roy-predicted real parts with the data (this 
corresponds to what we call consistency), they found that these constrain the 
two S-wave scattering lengths to the range Oq = 0.15 ±0.07, Oq = —0.053 ±0.028. 
In their subsequent work ||, they then used the Roy equations to solve the famous 
Up-Down ambiguity that occurs in the analysis of the S'-wave. 

The fact that, in their analysis, the matching point is taken below the mass of 
the p has an interesting mathematical consequence: As discussed in section 6.3, 
the Roy equations do then not admit a solution for arbitrary values of a{], Oq, even 
if cusps at the matching point are allowed for (the situation corresponds to row 
IV of table [l]). To enforce a solution, one may for instance keep the input for the 
imaginary parts as it is, but tune the scattering length a^. The result, however, 
in general contains strong cusps in the partial waves with I = 0, 1. These can 
only be removed if the input used for the imaginary parts above the matching 
point is also tuned - the situation is very different from the one encountered for 
our choice of the matching point. 

Basdevant, Froggatt and Petersen || |J constructed solutions of the Roy equa- 
tions by considering several phase shift analyses and a broad range of S'-wave 
scattering lengths. The method used by these authors is different from ours 
in that they relied on an analytic parametrization of the S- and P-waves from 
threshold up to \/s2 = a/110 M n = 1.47 GeV, the onset of the asymptotic region 
in their case. Some of the parameters occurring therein are determined from a fit 
to the data, some by minimizing the difference between the right and left hand 
sides of the Roy equations in the region below y/so = V&) M n = 1.08 GeV. In 
this manner, they construct universal bands corresponding to the Berkeley [|48H , 
Saclay ||5| and CERN-Munich phases as determined by Estabrooks et al. [|50]. 



The individual bands are not very much broader than the shaded region in fig. [Ti 
but they are quite different from one another: Crudely speaking, the Berkeley 
band is centered at the upper border of our universal band, while the one con- 
structed with the CERN-Munich phases is centered at the lower border. The 
Saclay band runs outside the region where we can find acceptable solutions at 
all. 
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In order to compare their results with ours, we first note that, for the six 
explicit solutions given in table 5 of @, the value of a[j varies between —0.06 and 
0.59. Only two of these correspond to values of the S'-wave scattering lengths in 
the region considered in the present paper: BKLY2 and SAC2. For these two, 
the value of the P-wave phase shift at 0.8 GeV is 108.3° and 108.0°, respectively, 
remarkably close to the central value of the range allowed by the data on the 
form factor, eq. (7.2). Concerning the value of 5® at 0.8 GeV, however, the two 
solutions differ significantly: While BKLY 2 yields 79.7° and is thus within our 
range in eq. (7.4), the value 70.2° that corresponds to SAC2 is significantly lower. 
In our opinion, that solution is not consistent with the experimental information 
available today. In the interval from threshold to 0.8 GeV, our solution differs 
very little from BKLY 2 . Above this energy, the imaginary part of the 1 = 
S'-wave in BKLY 2 is substantially smaller than the one we are using as an input. 
Nevertheless, the solutions are very similar at low energies, because the behaviour 
below the matching point is not sensitive to the input above 1 GeV. 

17 Summary and conclusions 

The Roy equations follow from general properties of the irn scattering amplitude. 
We have set up a framework to solve these equations numerically. In the following, 
we summarize the main features of our approach and the results obtained with it, 
omitting details - even if these would be necessary to make the various statements 
watertight. 

1. In our analysis, three energies So < s\ < s 2 play a special role: 

= 0.8 GeV , s = 32.9 Ml , 
y/s[ = 1.15 GeV , s 1 = 68 Ml , 
= 2 GeV , s 2 = 205.3 Ml . 

We refer to the point so as the matching point: At this energy, the region where we 
calculate the partial waves meets the one where we are relying on phenomenology. 
The point Si indicates the upper end of the interval on which the Roy equations 
are valid, while s 2 is the onset of the asymptotic region. 

2. Given the strong dominance of the S- and P-waves, we solve the Roy 
equations only for these, and only on the interval 4M^ < s < So, that is on the 
lower half of their range of validity. In that region, the contributions generated by 
inelastic channels are negligibly small. There, we set i]q(s) = r]l(s) = i]q(s) = 1. 
In the interval from sq to s 2 , we evaluate the imaginary parts with the available 
experimental information, whereas above s 2 , we invoke a theoretical represen- 
tation, based on Regge asymptotics. We demonstrate that crossing symmetry 
imposes a strong constraint on the asymptotic contributions, which reduces the 
corresponding uncertainties quite substantially - in most of our results, these are 
barely visible. 
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3. The Roy equations involve two subtraction constants, which may be iden- 
tified with the two S- wave scattering lengths a[j, Oq. In principle, one subtraction 
would suffice: The Olsson sum rule relates the combination 2 to an inte- 
gral over the imaginary parts in the forward direction (or, in view of the optical 
theorem, over the total cross section). This imposes a correlation between the in- 
put used for the imaginary parts and the values of the S'-wave scattering lengths, 
but using this constraint ab initio would lead to an unnecessary complication 
of our scheme. We instead treat the two subtraction constants as independent 
parameters. The consequences of the Olsson sum rule are discussed below. 

4. Unitarity converts the Roy equations for the S- and P-waves into a set 
of three coupled integral equations for the corresponding phase shifts: The real 
part of the partial wave amplitudes is given by a sum of known contributions 
(subtraction polynomial, integrals over the region sq < s < s 2 and driving terms) 
and certain integrals over their imaginary parts, extending from threshold to Sq. 
Since unitarity relates the real and imaginary parts in a nonlinear manner, these 
equations are inherently nonlinear and cannot be solved explicitly. 

5. Several mathematical properties of such integral equations are known, and 
are used as a test and a guide for our numerical work. In particular, the existence 
and uniqueness of the solution is guaranteed only if the matching point so is taken 
in the region between the place where the P-wave phase shift goes through 90° 
and the energy where the 1 = S'-wave does the same. As this range is quite 
narrow (0.78 GeV < y/so < 0.86 GeV), there is little freedom in the choice of the 
matching point - we use y/s^ = 0.8 GeV. According to table [I], the multiplicity 
index of the interval 0.86 < y/so < 1 GeV is equal to 1. By way of example 
{\/so = 0.88 GeV), we have verified that our framework indeed admits a one- 
parameter family of numerical solutions if the matching point is taken in that 
energy range. 

6. A second consequence of the mathematical structure of the Roy equations 
is that, for a given input and for a random choice of the two subtraction constants, 
the solution has a cusp at s^: In the vicinity of the matching point, the solution in 
general exhibits unphysical behaviour. The strength of the cusp is very sensitive 
to the value of a,Q. In fact, we find that the cusp disappears in the noise of 
our calculation if that value is tuned properly. Treating the imaginary parts as 
known, the requirement that the solution is free of cusps at the matching point 
determines the value of function of a[j. This is how the universal curve of 
Martin, Morgan and Shaw manifests itself in our approach. 

7. The input used for the imaginary parts above the matching point is subject 
to considerable uncertainties. In our framework, the values of the S- and P-wave 
phase shifts at the matching point represent the essential parameters in this 
regard. In order to pin these down, we first make use of the fact that the data on 
the pion form factor, obtained from the processes e + e~ — ► 7r + 7r~ and r — > tt~it i> t , 
very accurately determine the behaviour of the P-wave phase shift in the region 
of the p -resonance, thus constraining the value of 5;}(so) to a remarkably narrow 
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range. Next, we observe that the absolute phase of the irir scattering amplitude 
drops out in the difference 5\ — 5$, so that one of the sources of systematic 
uncertainty is eliminated. Indeed, the phase shifts extracted from the reaction 
ttN — > tcttN yield remarkably coherent values for this difference. Since the P- 
wave is known very accurately, this implies that Oq(so) is also known rather well. 
The experimental information concerning 8q, on the other hand, is comparatively 
meagre. We vary it in the broad range shown in fig. ||. 

8. The uncertainties in the experimental input for the imaginary parts and 
those in the driving terms turn the universal curve into a band in the (dp, dp) 
plane, part of which is shown in fig. [7[ Outside this "universal band", the Roy 
equations do not admit physically acceptable solutions that are consistent with 
what is known about the behaviour of the imaginary parts above the matching 
point. 

9. One of the striking features of the solutions is that, above the matching 
point, they very closely follow the real part of the partial wave used as input for 
the imaginary part, once the value of dp is in the proper range. The phenomenon 
is discussed in detail in section [H], where we show that, in a certain sense, this 
property represents a necessary condition for the solution to be acceptable phys- 
ically. The region where this consistency condition holds is shown in fig. |9]: It 
roughly constrains the admissible values of a% to the lower half of the universal 
band. 

10. As mentioned above, the Olsson sum rule relates the combination 2d[j— 5dg 
of scattering lengths to an integral over the imaginary parts of the amplitude. 
Evaluating the integral, we find that the sum rule is satisfied in the band spanned 
by the two red curves shown in fig. ^. The Olsson sum rule thus amounts to 
essentially the same constraint as the consistency condition. Presumably, the 
universal band is of the same origin: Physically acceptable solutions only exist 
if the subtraction constants are properly correlated with the imaginary parts. 
The shaded region in fig. |9] shows the domain where all of these conditions are 
satisfied. It is by no means built in from the start that the various requirements 
can simultaneously be met - in our opinion, the fact that this is the case represents 
a rather thorough check of our analysis. 

11. The admissible region can be constrained further if use is made of ex- 
perimental data below the matching point. At the moment there are two main 
sources of information on tttc scattering below 0.8 GeV: A few data points for the 
1 = 2 S'-wave phase shift - which to our knowledge will, unfortunately, not be 
improved in the foreseeable future - and a few data points on 5® — 5\ very close 
to threshold, as measured in K e ^ decays. These data do provide an important 
constraint. We compare our solutions inside the universal band to both sets of 



data. As shown in fig. 12, the corresponding \ contours nicely fit inside the 



universal band. The net result for the allowed range of the parameters is shown 
in fig. [L7], which summarizes our findings. 

12. To our knowledge, the Roy equation analysis is the only method that 
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allows one to reliably translate low energy data on the scattering amplitude into 
values for the scattering lengths. As discussed above, the available data do cor- 
relate the value of Oq with the one of a{]. Unfortunately, however, the value of a[j 
as such is not strongly constrained: In agreement with earlier analyses, we find 
that these data are consistent with any value of a{] in the range from 0.18 to 0.3. 

13. The new experiments at Brookhaven p7| and at DA$NE will yield 
more precise information in the very near future. We expect that the analysis of 
the forthcoming results along the lines described in the present paper will reduce 
the error in Oq by about a factor of three. Moreover, the pionic atom experiment 
under way at CERN JZD) will allow a direct measurement of — a^\ and thus 
confine the region to the intersection of the corresponding, approximately vertical 
strip with the region shown in fig. [T7. 

14. The two subtraction constants the essential parameters at low 
energies: If these were known, our method would allow us to calculate the S- 
and P-wave phase shifts below 0.8 GeV to an amazing degree of accuracy. The 
parameters Oq, ag act like a filter: If the solutions are sorted out according to 
the values of these parameters, the noise due to the uncertainties in our input 
practically disappears, because variations of that input require a corresponding 
variation, either in a[] or in - otherwise, the behaviour of the solution near 
the matching point is unacceptable. A simple explicit representation for the S- 
and P-wave phase shifts as functions of the energy is given in appendix [D]. The 
representation explicitly displays the dependence on a(j, a^. 

15. We have also analyzed the implications for the scattering lengths of the 
P-, D- and P-waves, as well as for the various effective ranges. The fact that 
aj] and <2q are the essential low energy parameters manifests itself also here: If 
we change the input in the Roy equations within the uncertainties, but keep Oq 
and dp constant, the values of the various threshold parameters vary only by tiny 
amounts, typically around one percent or less. The main source of uncertainty in 
the determination of the threshold parameters is by far the one attached to the 
.S-wave scattering lengths. 

16. If the energy approaches the matching point, the uncertainties in the 
experimental input, naturally, come more directly into play. Also, the uncertain- 
ties in the driving terms grow rather rapidly with the energy. At the kaon mass, 
however, these are still very small. We have analyzed the phase shifts at E = Mk 
in detail, because these represent an important ingredient in the calculation for 
various decay modes of the K mesons. The result shows that the uncertainties 
are dominated by those in a^, Oq, also at that energy. We conclude that the fu- 
ture precision data on P^ 4 -decay and on pionic atoms will translate, via the Roy 
equations, into a rather precise knowledge of the tttc scattering amplitude (not 
only the lowest three partial waves) in the entire low-energy region, extending 
quite far above threshold. 

17. In the present paper, we followed the phenomenological path and avoided 
making use of chiral symmetry, in order not to bias the data analysis with the- 
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oretical prejudice. A famous low energy theorem |32] predicts the values of the 
two basic low energy parameters in terms of the pion decay constant. The predic- 
tion holds to leading order in an expansion in powers of the quark masses. The 
corrections arising from the higher order terms in the chiral expansion are now 
known to order p 6 (two loops) [[H]]. We plan to match the chiral perturbation 
theory representation of the scattering amplitude with the phenomenological one 
obtained in the present paper f80fl . This should lead to a very sharp prediction 



for <2q and a 2 ,. The confrontation of the prediction with the forthcoming results 
of the precision measurements will subject chiral perturbation theory to a crucial 
test. 
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A Integral kernels 

Crossing symmetry, A(s, u,t) = A(s,t,u), implies that the isospin components 
T = (T°, T 1 , T 2 ) are subject to the constraints [u = 4M 2 — s — t) 

T(s,u) = C tu f(s,t) , 
f{t,s) = C st f{s,t), 
T(u,t) = C su T(s,t) , 

where the crossing matrices C tu = C ut , C su = C us , C st = C ts are given by 
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Their products obey the relations 

{Ctu) 2 — (C su ) 2 = (C st ) 2 — 1 ) 

St ^tU — ^tU ^US — ^US ^St ; <-> Ut ^ SU ~ <- ' Ut ^ts ■ 
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The quantities #2(5, t, s'), ^(s, t, s') occurring in the fixed-t dispersion relation 
represent 3x3 matrices built with C st , C tu and C su , 

i /in 
s,t,s) = — — -—p- [uC st + sC st C tu ) 



AMI) Vs'-t s'-« 

M') = T^-A^— + ^ L -) , (A.l) 

7TS'{S —Uo) \S ' — S S' — UJ 

where u = AMI - s-t and u = AMI ~ t- 

The straightforward partial wave projection of the amplitude reads 

t I t (s) = -^j\zP e (z)T I (s,t z ) , t z = \(AMl-s){l-z) . (A.2) 

On account of crossing symmetry, the formula is equivalent to 

4(s) = ^r ( X dzP l {z)T'{a t U) . (A.3) 

3Z7T JO 

As pointed out by Roy the second form of the projection is preferable in the 
present context, because it involves smaller values of \t z \, so that the domain of 
convergence of the partial wave series for the imaginary parts on the r.h.s. of the 
fixed-t dispersion relation (|2.4j) becomes larger: Whereas for the straightforward 
projection, the large Lehmann-Martin ellipse is mapped into — AM% < s < 32M%, 
the one in eq. ( |A.3| ) corresponds to — AM% < s < 60M^. 



The kernels K^,'(s,s') that occur in eq. (|1.1| ) are different from zero only if 
both I+£ and I'+£' are even. With the partial wave projection ( |A.3| ), the explicit 
expression becomes^ 

K I J(s,s') = (2£'+l) f 1 dzPe(z)K e ,(s,t z ,s') n ' , 

J 

t z = \{AMl-s){l-z) . (A.4) 
The functions Kp(s, t, u) 11 ' are the matrix elements of 

K e ,{s, t, s') = g 2 (s, t, s') + g 3 (s, t, s') pJl+ - _ 2 ^ 2 J . (A.5) 



llT 



Note that the fixed-t dispersion relation (2.4) is not manifestly crossing symmetric - for 



£' > 2, the kernels do depend on the specific form used for the partial wave projection. In 
particular, the kernels occurring in the Roy equations for the waves with £ > 3 are proportional 



to (s — 4M%Y only if the projection in eq. (A.2) is used - for the one we are using here, the 



proper behaviour of the solutions only results if the contributions from the imaginary parts 



of the different partial waves compensate one another near threshold (see section 14.2). For a 
detailed discussion of these issues we refer to fn\. 
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The kernels contain the usual pole at s = s', generating the right hand cut 
of the partial wave amplitudes, as well as a piece K\\, (s, s') that is analytic in 
Res > 0, but contains a logarithmic branch cut for s < — (s' — 4M 2 ): 

Kg {s, s>) = ^J—^ 5 11 ' 5 U , + Kg (s, s') . 

To illustrate the structure of the second term, we give the explicit expression for 
I = I' = £ = £' = 0: 

K™(s s>) - 2 fJ S + S '~ AM2 A - 2s + 5s '- 16M - 

oo[S,S) ~ 3n(s-4Mi) m { s' J ^s'{s'-AMl) ' 

We do not need to list other components - they may be generated from the above 
formulae by means of standard integration routines. 



B Background amplitude 

B.l Expansion of the background for small momenta 

The background amplitude only contains very weak singularities at low energies. 
At small values of the arguments, A(s, t, u)d thus represents a slowly varying 
function of s, t, u, which is adequately approximated by a polynomial. We may, 
for instance, consider the Taylor series expansion around the center of the Man- 
delstam triangle: Set s = |M^, s = so + x, t = s — |(x — y), expand in powers 
of x and y and truncate the series. Alternatively, we may exploit the fact that, 
in view of the angular momentum barrier, the dispersion integral over the imag- 
inary parts of the higher partial waves receives significant contributions only for 
s'>lGeV 2 . For small values of s and t, we can therefore expand the kernels 
g2{s, t, s') and gs(s, t, s') in inverse powers of s'. The coefficients of this expansion 
are homogeneous polynomials of s, t and M 2 , which may be ordered with the 
standard chiral power counting. The corresponding expansion of the Legendre 
polynomial starts with 

Truncating the expansion at order p 6 , the background amplitude becomes 

T(s,t) d = -32n {(tuC st + suC su + stC tu ) (1 + C SU ) T (B.l) 

+{s 2 1 C tu + u 2 s C su + t 2 u C st + (t 2 s C tu + s 2 u C su + u 2 1 C st ) C su } h 
+stu(l + C su )H} + 0(p 8 ) . 
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The coefficients Iq and I\ represent moments^ of the imaginary part at t = 0, 

ri ^ 1 r ds Imr f (s,0) d 

n 32vr 2 i4Af2 S «+2( S -4M2) ' 1 ' ' 

In view of the optical theorem, these quantities are given by integrals over the 
total cross section, except that the contributions from the S- and P-waves below 
s 2 are to be removed. Equivalently, we may express these coefficients in terms of 
the imaginary parts of the partial waves: 

j " (2i + l) /"' dslmt^s) " (2i + l) /-°° dslmti(s) 

n f^ 2 it Jam? s n + 2 (s - 4M2) + ^ tt J S2 s n + 2 (s - 4M 2 ) ' 1 ' ' 

Except for a contribution proportional to l\, the last term in eq. QB.1|) may be 
expressed in terms of the derivative of Im T(s,t)a with respect to t: 

1 1 32 vr 2 hul s 3 dt t =o 

Here, only the higher partial waves contribute: 

H' = ±(21 + l) W + 1) - 2 1 r sr ^l 1 ^ • (B-5) 

The expression is similar to the one for l[, except that the sum over the angular 
momenta picks up a factor of t{t + 1), indicating that partial waves with higher 
values of I are more significant here. Note that all of the above moments are 
positive. 



B.2 Constraints due to crossing symmetry 

The expansion of the background amplitude starts at order p 4 , with a manifestly 
crossing symmetric contribution determined by the moments Iq. The term from I\ 
is also crossing symmetric, but the one proportional to s t u violates the condition 
T(s,u)d = C tu T(s,t)d, unless the 1 = 1 component of the vector (1 + C su ) H 
vanishes, i.e. 

2H° = 9H 1 + 5H 2 . (B.6) 

This sum rule is both necessary and sufficient for the polynomial approxima- 
tion to the background amplitude to be crossing symmetric up to and including 
contributions of order p 6 . 

12 The factor l/(s — 4M 2 ) could also be expanded in inverse powers of s, but this would 
worsen the accuracy of the polynomial representation. Note that the same factor also occurs 
in the representation ([T^) for the contributions generated by the imaginary part of the S- and 
P-waves below S2'- The expansion of the functions W 1 (s) in powers of s yields integrals of the 
same form. Hence the low energy expansion of the full amplitude can be expressed in terms of 
moments of this type. 
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The sum rule illustrates the well-known fact that crossing symmetry leads 
to stringent constraints on the imaginary parts of the partial waves with t > 2 
(for a thorough discussion, see |H, |4!|). Crossing symmetry implies for instance 



that Imt^s) can be different from zero only if some of the higher partial waves 
with I — 1 or / = 2 also possess an imaginary part - in marked contrast to the 
situation for the S- and P-waves, where crossing symmetry does not constrain 
the imaginary parts. 

In the form given, the sum rule only holds up to corrections of order M%. 
We may, however, establish an exact variant by expanding the 1 = 1 component 
of the relation T(s,u)d = C 4u T(s,t)d around threshold, for instance in powers 
of t and u. In order for the term of order tu occurring in the expansion of the 
left hand side to agree with the corresponding term on the right hand side, the 
imaginary parts must obey the sum rule 

rOO fig , 

L s^s-mi) {21-^(..0)-5ImI«(.,0)} = 

3 £; {(' - 4Ag)Imf (..0) - 2Imr-(,,0)} , (B.7) 

where T 7 (s,t) stands for the partial derivative of T 7 (s,t) with respect to t. Ex- 
pressed in terms of the imaginary parts of the partial waves, the relation reads 

/•OO r] a 

E <*+ D jL, ?(s^Ml? {2Imi?(s) - 5Im « s » = 

£ (M + 1) {t (t + l) - 2} £ g/;;^ 91-^) . (B.8) 



£=2,4, 



£=3,5 



The approximate version ( [B.5| ) differs from this exact result only through terms 
of order M%. 

The constraints imposed by crossing symmetry show, in particular, that the 
concept of tensor meson dominance is subject to a limitation that does not occur 
in the case of vector dominance: The hypothesis that convergent dispersion inte- 
grals or sum rules are saturated by the contributions from a spin 2 resonance leads 
to coherent results only at leading order of the low energy expansion. The sum 
rule ( |B.7| ) demonstrates that the hypothesis in general fails: Crossing symmetry 
implies that singularities with i > 2 cannot be dealt with one by one. 

Since the relation ( |B.6|) ensures crossing symmetry, the above low energy 
expansion of the isospin components of the amplitude is equivalent to a manifestly 
crossing symmetric representation of the background amplitude: 

A(s, t, u) d = pi + p 2 s + p 3 s 2 + p 4 (t - uf +p 5 s 3 + p 6 s(t - uf + 0(p 8 ) . (B.9) 

By construction, A(s,t,u)d does not contribute to the S"-wave scattering lengths. 
This condition fixes pi and p 2 in terms of the remaining coefficients: 

P! = -16M> 4 , p 2 = 4M^(-p 3 +p4-4M^p 5 ), (B.10) 
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The explicit expressions for the latter read 

ps = ^ « - Wo 1 - i%) + ^fK (-8 A° - 21A 1 + ml + 12 H) , 

P4 = S^{ll + ll) + lQixMl{ll + ll) , (B.ll) 
p 5 = ^(8/ 1 ° + 9/ 1 1 -ll/ 1 2 -6iJ) , 
p 6 = 47r(/ 1 1 -3J 1 2 + 2#) . 

In view of the sum rule (p.6j ), only two of the components of H are independent. 
Moreover, the amplitude only involves a combination thereof: 

H = §(#° - 2H 1 ) = I(H° + 2H 2 ) = H 1 + H 2 . (B.12) 

The above formulae show that the leading background contribution is deter- 
mined by the integrals Iq , which yield 

Pl = 0(Mi) , V 2 = O(Ml) , p 3 = 0(1) , p 4 = 0(1) . 

The contributions from Ii and H modify the result by corrections that are sup- 
pressed by one power of M 2 and, in addition, generate a polynomial of third 
degree in s, t, u, characterized by p§ and pe- 

B.3 Background generated by the higher partial waves 

Next, we turn to the numerical evaluation of the integrals Iq, Ii, H and first 
consider the contributions from the imaginary parts of the partial waves with 
£ > 2 in the region below 2 GeV. The integrals are dominated by the resonances, 
which generate peaks in the imaginary parts. In the vicinity of the peak, we may 
describe the phase shift with the Breit-Wigner formula 

2i5 { s) = M r 2 + iT r M r - s 
' M? - iT r M r - s ' 

where M r and T r denote the mass and the width of the resonance, respectively. 
To account for inelasticity (decays into states other than tttt), we multiply the 
corresponding expression for the imaginary part of the partial wave amplitude 
with the branching fraction IV^^/TV. This leads to 

s r r M 2 

S ~ AMI - M r ) 2 + T r M r ' 

where I r and l r denote the isospin and the spin of the resonance, respectively. In 
the narrow width approximation, the formula simplifies to 

Imt£(a) = 7rr™M r (l - 4M 2 /M 2 )^6(s - M 2 ) . (B.13) 
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Only four of the states listed in the particle data booklet |70j below 2 GeV 
have spin I > 2 and carry the proper quantum numbers to be produced in tttt 
collisions: The spin 2 resonances /2(1275) and /2(1525), the spin 3 state ps(1681) 
and the state /j(1710), whose spin is not firmly established, but must be even. 
There is no evidence for exotic states: f 2 and fj are isoscalars, while the ps is 
an isovector. 

Very likely, the lightest spin 4 state is the / 4 (2044): A linear p(770) — 
/ 2 (1275) — p 3 (1691) Regge trajectory calls for a spin 4 recurrence almost ex- 
actly there. At any rate, if the spin of the state /j(1710) were equal to 4 or 
even larger, it would sit above that trajectory and thus upset the standard Regge 
picture, which we will be making use of to estimate the asymptotic part of the 
driving terms. We take it for granted that J = or 2 and conclude that only the 
1 = D-w&ve and the F-wave contain resonances below 2 GeV. In the follow- 
ing, we discuss the contributions generated by these states, comparing the result 
obtained from the narrow width formula with the one found on the basis of two 
different phase shift analyses. 

The most important contribution arises from the tensor meson / 2 (1275). In- 
serting the values Mf 2 = 1275 MeV, Tf 2 ^ wn = 157 MeV, the narrow width formula 
gives J ° /2 = .25 GeV~ 4 , l[ h = .15 GeV -6 , H% = .93 GeV~ 6 , to be compared with 
the results obtained with the parametrizations of the D-wave in refs. [[47}] an d [ [57H , 
which yield 

[47] : J ° D = .25 GeV" 4 , 1% = .18 GeV" 6 , H° D = 1.10 GeV" 6 , (B.14) 
: J ° D = .27GeV- 4 , 1% = .19GeV" 6 , H° D = 1.17 GeV" 6 . (B.15) 

These numbers show that the contributions from the imaginary part of the D- 
wave are dominated by the /2(1275). 

We add a few remarks concerning the detailed behaviour of Imt^s) and first 
note that the /2(1525) mainly decays into KK. In the present context, this state 
may be ignored, because the corresponding tttt partial width is tiny: Tf^^ = 
.62±.14 MeV. The phase shift analysis of ref. [|57]] does contain a second resonance 
in the D-wave, which generates a small enhancement in the integrands on the 
r.h.s. of eqs. (|B.3| ), ( p.5| ) towards the upper end of the range of integration. The 
numerical result in eq. ( |B.15[ ) includes the tiny contribution produced by this 
enhancement, but this effect only accounts for a small fraction of the difference in 
the values obtained with the two different phase shift analyses. The main reason 



for that difference is that the two representations of the D-wave in refs. p7| , [57 
do not agree very well on the left wing of the /2(1275). In the context of the 
present paper, these details are not essential - we use the difference between the 
two phase shift analysis as a measure for the uncertainties to be attached to the 
moments. 

To estimate the significance of the remaining partial waves with / = 0, we 
consider the contribution generated by the /4(2044). This resonance also mostly 
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decays into states other than tttt. The relevant partial width is r^-^ = 35 ± 
4MeV. The narrow width formula shows that the contribution from this state is 
very small: J ° /4 = .009GeV~ 4 , I% u = .002GeV~ 6 , H% = .04GeV~ 6 . Moreover, 
the center of the peak is outside our range of integration - more than half of the 
contribution from this level is to be booked in the asymptotic part. We conclude 
that the imaginary parts of the partial waves with i > 4 only matter at energies 
above 2 GeV. 

The p3(1681) shows up as a peak in the imaginary part of the F-wave. Accord- 
ing to the particle data tables [ff5| , it mainly decays into 47r. The partial width of 
interest in our context is r p3 ^ n7T = 38 ±3 MeV. Inserting this in the narrow width 
formula, we obtain I^ p3 = .020 GeV" 4 , l\ p:i = .007 GeV" 6 , Hj 3 = .07 GeV" 6 , to 
be compared with the values found by performing the numerical integration with 
the representations for the F-wave given in the two references quoted above: 

[47] : ll F = .028 GeV~ 4 , l{ F = .012 GeV" 6 , H l F = .12 GeV" 6 , (B.16) 
|57j : I] F = .030 GeV~ 4 , l\ F = .016 GeV" 6 , H F = .16 GeV" 6 . (B.17) 



In the present case, the narrow width formula only accounts for about half of the 
result: The region below the resonance is equally important. There, the difference 
between the two phase shift analyses is more pronounced than for the F-waves. 
Accordingly, the uncertainties in the F-wave contributions to the moments are 
larger. 

The formula ( [B.13| ) predicts that the contribution generated by the imagi- 



nary part of the 1 = 2 waves vanishes, because that channel does not contain 



any resonances. According to Martin, Morgan and Shaw |82|| , the F-wave phase 
shift may be approximated as 5f(s) ~ —0.003 (s/4M^) (1 — 4M^/s)a. The corre- 

T2 



sponding contributions to the moments are indeed very small: Iq = 0.005 GeV 4 , 



If = 0.006 GeV -6 , H = 0.04 GeV" 6 . In the following, we assume that these 
estimates do hold to within a factor of two. 

This completes our discussion of the contributions generated by the higher 
partial waves in the region below 2 GeV. The net result is that these are due 
almost exclusively to the D- and F-waves. The numerical results are listed in 
row L of table |6]. For I = 0, 1, the values given rely on the phase shift analyses 
of refs. p7| , |57f , while the estimates for 1 = 2 correspond to the parametrization 
of ref. 
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/ = 


I = 1 


1 = 2 








H° 




n 


H 1 


T 2 


n 


H 2 




GeV -4 


GeV -6 


GeV -6 


GeV -4 


GeV -6 


GeV -6 


GeV -4 


GeV -6 


GeV -6 


L 


.26 


.19 


1.13 


.029 


.014 


.14 


.005 


.006 


.04 


R 


.03 


.004 


.11 


.018 


.003 


.07 








P 


.01 


.001 


.04 


.010 


.001 


.04 


.010 


.001 


.04 


total 


.30 


.19 


1.28 


.058 


.018 


.24 


.015 


.007 


.08 


± 


.01 


.01 


.05 


.007 


.002 


.03 


.008 


.006 


.04 



Table 6: Moments of the background amplitude. The rows L, R and P indicate 
the contributions from the region below 2 GeV, from the leading Regge trajectory 
and from the Pomeron, respectively. The last two rows show the result for the 
sum of these contributions and our estimate of the uncertainties, respectively. 

B.4 Asymptotic contributions 

We now turn to the contributions from the high energy tail of the dispersion 
integrals. The asymptotic behaviour of the scattering amplitude may be analyzed 
in terms of Regge poles. A trajectory with isospin I generates a contribution 
oc s a ^ to the t-channel isospin component ImT^(s,t), which is defined by 

ImT(')(M) = E^/' ImT/ 'M) • 
v 

The asymptotic behaviour of the amplitude with I t — 1 (s — > oo, t fixed) is 
governed by the p -trajectory, 

ImT (1) (s,t) =f3 p {t)s a ^ t) . 

The Pomeron dominates the high energy behaviour of the I t = amplitude. 
Together with the contribution from the /-trajectory, the Regge representation 
of this component reads 

Im T (0) (s, t) = 3 P(s,t) + P f (t) s a ' w . 

In the absence of exotic trajectories, the component with I t — 2 rapidly tends to 
zero when s becomes large. The asymptotic behaviour of the s-channel isospin 
components thus takes the form 

Im T°(s, t) = P(s, t) + \(3 f {t) s a ' (t) + (3 p (t) s ap{t) + (t <-> u) , 
lmT l (s,t)=P(s,t) + l(3 f (t)s a ^ + \(5 p {t)s a "^ - (t ^ u) , (B.18) 
Im T 2 (s, t) = P{s, t) + \(3 f (t) s a f {t) - \(5 p {t) s a » (t) + (t <-> u) . 

If t is kept fixed, the terms with P(s, t) and /3(t) s a ^ dominate, generating a peak 
in the forward direction, while the analogous structure in the backward direction 
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(fixed u) is described by those with P(s, u) and /3(u) s a ( u >. At fixed t, the crossed 
terms drop off very rapidly with s, so that their contribution disappears in the 
noise of the calculation and may just as well be dropped. 

The Lovelace-Shapiro- Veneziano model |83|, |84], |85|] provides a very instructive 
framework for understanding the interplay of the asymptotic contributions with 
the resonance structures seen at low energies (see appendix H). In that model, 
the p- and /-trajectories are linear and exchange degenerate, 

ot p (t) = otf(t) = a + t q;i . (B.19) 

We fix the intercept with the Adler zero, at(M%) = |, and choose the slope such 
that the spin 1 state on the leading trajectory occurs at the proper mass: 

ttl = \ (M 2 p - Ml)- 1 , a = \ - ai Ml . (B.20) 

The amplitude may be represented as a sum of narrow resonance contributions. 
Since the model does not contain exotic states, ImT 2 (s,t) vanishes, so that the 
residues and f3 p {t) are in the ratio 3:2. The explicit expression reads 

««> = = 4rff ■ (R21) 

Finally, we fix the overall normalization constant A such that the width of the p 
agrees with what is observed. This requires 

A = 96 7r T p M 2 p (Mp — 4Ml)~^ . (B.22) 

The model explicitly obeys crossing symmetry and yields a decent picture 
both for the masses and widths of the resonances occurring on the leading tra- 
jectory and for the qualitative properties of the Regge residues (3 p (t), (3 fit). The 
main deficiency of the model is lack of unitarity: It does not contain a Pomeron 
term, so that the total cross section tends to zero at high energies. While the 
model yields quite decent values for the full widths, it does not account for the 
fact that the resonances often decay into states other than tttt, particularly if 
the available phase space becomes large - in the model, the branching fraction 
r r _> n7T /r r is equal to 1. Consequently, the LSV-model overestimates the magni- 
tude of the Regge residues - a significant fraction thereof should be transferred 
to the Pomeron term. For this reason, the model can only serve as a semi- 
quantitative guide. 

As discussed in section B.2| , crossing symmetry strongly correlates the asymp- 
totic behaviour of the partial waves with their properties at low energy. In par- 
ticular, the parameters occurring in the Regge representation of the scattering 
amplitude can be extracted from low energy phenomenology. For a review of 
these calculations, we refer to the article by Pennington |[43||. The value obtained 
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for /3 P (0) is smaller^ than what follows from eqs. ( |B.21| ), ( |B.22| ) by a factor of 



0.6 ±0.1. Also, while the formula (|B.21|) implies that the residue contains a zero 
at t = 2M 2 — Mp = —0.55 GeV 2 because a(t) vanishes there, the calculation 
of ref. [43| instead yields a zero at t = -0.44 ± 0.05 GeV 2 . This confirms the 



remarks made above: The LSV-model describes the qualitative properties of the 
Regge residues quite decently, but overestimates their magnitude. 

In the numerical evaluation, we use the linear p -trajectory specified above, 
a p (t) = a(t), and fix the corresponding residue with the results of ref. which 
are adequately described by a modified version of the LSV-formula: 

W = rll ia + \ 1 > to = -°' 44 GeV2 ' Al = (' 78 ± - 13 ) A • ( R23 ) 
T[(t -tojaij 

We determine the properties of the /-trajectory with exchange degeneracy, i.e. set 
Oif{t) = a(t) and /3f(t) = |/3 p (t). For the Pomeron, we use the representation 

P{s,t) =ase^ bt . (B.24) 



While the parameter 6 = 8 GeV -2 j|3| describes the width of the diffraction 
peak, the optical theorem implies that a represents the asymptotic value of the 
total 7T7T cross section. Evidently, the above parametrization can be adequate 
only in a limited range of energies: The cross section does not tend to a constant, 
but grows logarithmically. In the present context, however, the behaviour at very 
high energies is an academic issue, because the integrands of the moments rapidly 
fall off with s. What counts is that the above representation yields a decent 
approximation for cm. energies in the range between 2 and 3 GeV. There, the 
terms generated by the p -/-trajectory are by no means negligible: The formula 
(P.18| ) shows that at 2 GeV (3 GeV), these terms by themselves generate a 
contribution to ImT°(s, 0) that corresponds to a total cross section of 21 mb (14 
mb) - in the energy range relevant for the moments, the Pomeron term does not 
represent the dominating contribution to the total cross section. As discussed 



in detail in ref. [ 43] , crossing symmetry leads to the estimate o = (6 ± 5) mb. 



Although the error bar is large, the value is significantly smaller than what is 
indicated by the rule of thumb a™ t ~ | ~ | a^J 1 ~ 20 mb. 

Indeed, the sum rule ( |B.6| ) confirms this result. The numerical values obtained 
with the above representation for the contributions from the p -/-trajectory are 
indicated in row R of table If the high energy tail is omitted altogether, the 
l.h.s. of the sum rule ( [B.6] ) becomes (2H°)l = 2.3 GeV -6 , while the r.h.s. amounts 
to (9H 1 + 5H 2 ) L = 1.5 GeV -6 . Clearly, further contributions are required to 
bring the sum rule into equilibrium. The Regge terms do contribute more to the 



13 In ref. Ji]|, the residue is written as (3 p (t) = ^Tr^ p (t) a" p ' 2 . The result obtained for the 
value at t = is 7 P (0) = (0. 6 ± O.l^ M" 1 , to be compared with the number 7 P (0) = 0.97 M" 1 



that follows from eqs. ( |B.19D -( |B.22D 
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right than to the left and reduce the discrepancy by a factor of two. Since the 
Pomeron affects the various isospin components almost equally, it contributes 
about 7 times more to the right than to the left. For the sum rule to be obeyed 
within the uncertainties of the remaining contributions, the value of a must be 
in the range a = (5 ± 3) mb. 

Let us compare our representation of the background with the model used 
for the asymptotic behaviour in the early literature. Assume that, above an 
energy of 1.5 GeV, the imaginary parts can be described by a Pomeron term with 
Otot = 20 mb and a Regge term that corresponds to the leading trajectory of the 
LSV-model. The l.h.s. of the sum rule (|B.6|) then takes the value 2H° = 3.3, while 
the r.h.s yields 9H 1 + 5H 2 = 6.1 (to be compared with the value 2.6 obtained for 
either one of the two sides with our representation of the background). Evidently, 
the model is in conflict with crossing symmetry. In the region relevant for the 
driving term integrals, the LSV-model overestimates the Regge residues by about 
40% j|3| and the sum rule ( |B.6|) then implies that the value a = 20 mb is too 



large by about a factor of 4. 

We repeat that our calculation has no bearing on the asymptotic behaviour 
of the total cross section - we are merely observing that, unless the value of a is 
in the range 5 ± 3 mb, the representation used for the amplitude violates crossing 
symmetry. The row P indicates the contributions to the moments generated by 
the Pomeron if o is taken in the middle of this range. The net result of our 
calculation is contained in the last two rows of table || which list the outcome 
for the moments and for the error bars to be attached to these, respectively. For 
the quantity H defined in eq. ( p,12| ), we obtain 



H = 0.32 ± 0.02 GeV" 6 . (B.25) 



B.5 Driving terms 

The polynomial approximation for the background amplitude can be used to 
determine the low energy behaviour of the driving terms - it suffices to evaluate 
the partial wave projections of the polynomial T(s, t)d- The range of validity 
of the resulting representation for the driving terms, however, only extends to 
cm. energies of about 0.6 GeV. For our numerical work, we need a representation 
that holds for higher energies. 

The approximations for the imaginary parts discussed above yield the follow- 
ing representation of the driving terms: 

d\{s) = 4(s) L + d{(s) R + d\(s) P , 
1 f 1 

dI e(s) H = — / dzP^T^s^H , H=R,P 

62, 7T JO 
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T(s,t) H = ds'g 2 (s,t, s') ■lmT(s', 0) H + / ds'g 3 (s, t, s') ■ Im T(s', t) H , 

Im T°(s, t) fi = | (3 p (t) s«W + f /3» s«W , 
ImT 1 ^,^ = s«W - /3 p ( u ) S °W , 
lmT 2 (s,t) R = , 

ImT°(s, t) P = ImT 2 (s, t) P = P(s, t) + P(s, u) , 
lmT\s, t) P = P(s, t) - P(s, u) . 

The result of the numerical evaluation of these integrals with the parameter values 
specified above is given in eq. ( f4.1|) . 

We use the difference between the results for g?q(s)l and d\(s)i obtained with 
the two phase shift analyses quoted above as a measure for the uncertainties in 
these quantities. In the case of the 1 = 2 D-w&ve, we assume that the Martin- 
Morgan-Shaw formula does describe the behaviour of the imaginary part to within 
a factor of 2. For the Regge-contributions, we use the error estimate 7 P (0) = (0.6± 
O.^M" 1 given in ref. Finally, the uncertainties attached to the Pomeron 



term correspond to those in the value a = 5 ± 3 mb, obtained in section |B.4 
The result quoted in eq. (|4.2j ) is obtained by adding the corresponding error bars 
quadratically and fitting the outcome with a polynomial. 

There is a neat and rather thorough check of the above calculation. The 
driving terms represent the partial wave projections of the background amplitude. 
Since that amplitude must be crossing symmetric, we may equally well calculate 
the projections with the formula ( |A.2| ) instead of using ( |A.3| ) - the result should be 



the same. The modification of the partial wave projection changes the form of the 
kernels Kjj, (s, s') and the contributions from the imaginary parts of the higher 
partial waves below 2 GeV then change, quite substantially. The contributions 
from the asymptotic region, however, are also modified. In the sum, these changes 
indeed cancel out, to a remarkable degree of accuracy. This corroborates the 
claim that our description of the background is approximately crossing symmetric. 
Evidently, the sum rule ( |B.6|) plays an important role here, as it correlates the 



magnitude of the asymptotic contributions with those from the low energy region. 



C Sum rules and asymptotic behaviour 
C.l Sum rules for the P-wave 



As discussed in section [II], the Olsson sum rule ensures the correct asymptotic 
behaviour of the t-channel 1 = 1 scattering amplitude T^(s, t) for s — ► oo, t = 0. 
The requirement that this amplitude has the proper high energy behaviour also 
for t < implies a further constraint, which is readily derived from the fixed-t 
dispersion relation (|2.4j ). It suffices to evaluate the coefficient of the term that 
grows linearly with s and to subtract the value at t = 0. The result involves the 
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following integrals over the imaginary parts of the amplitude (t < 0): 

cm- r a 2ImT°(M) + 3Imf 1 M)-5ImT 2 M) fn 
S{t) = L/ S 12s(s + t-4M*) 

°° (s-2M2)ImT 1 (s,0) 



4M2 s(s-4M2)(s-t)(s + t-4M2) ' 

The barred quantities stand for 1m T I (s,t) = {ImT J (s,t) — ImT^s, 0)}/i. The 
amplitude T"W(s, £) has the proper asymptotic behaviour only if S(t) vanishes for 
space-like values of t. Since the 5-waves drop out, the condition amounts to a 
family of sum rules that relate integrals over the imaginary part of the P-wave 
to the higher partial waves. For t — 0, for instance, the sum rule may be written 
in the form 



2 



The integrals over the individual partial waves converge more rapidly than in the 
case of the Olsson sum rule, but the factor + 1) gives the higher partial waves 
more weight - in fact, the contributions from the asymptotic region are even more 
important here. The sum rule is of the same structure as the one that follows 
from crossing symmetry, eq. (B.7), but there are two differences: The integrals 
converge less rapidly by one power of s and the P-wave does not drop out. 

Since the sum rule ( |C2| ) offers a good opportunity to check the representation 



used for the asymptotic region, we evaluate it explicitly with our input for the 
imaginary parts. We split the integration into one from threshold to — 2 GeV 
and one over the asymptotic region, s > S2 (compare appendix B). Denoting the 
low energy part of the integral over the P-wave by 

s P =! S2 ds Imti(s) 



/4A/2 s 2 (s-AMl) ' 

we write the sum rule in the form 

S P = S D + S F + S as , (C.3) 

where S D and S F stand for the analogous integrals over the D- and P-waves. 
While the low energy contributions from the waves with I > 4 are neglected, their 
high energy tails are accounted for in the term S as , which collects all contributions 
from the region above S2- 

The form (|C.2|) of the sum rule is suitable to calculate the contributions from 
the interval 4M 2 < s < S2- Numerically, we obtain 

S P = 1.93 ±0.08 , S D = 0.55 ± 0.03 , S F = 0.13 ±0.01 , 



72 



in units of GeV 4 . To evaluate the asymptotic contributions, we instead use the 



form (|C.1|) : The term S as coincides with the expression 5(0)/487r, except that 
the integration now only extends over the interval S2 < s < 00. Inserting the 
representation specified in appendix B.4, we find that the bulk stems from the 
leading Regge trajectory (1.12 ± 0.19). The Pomeron does not contribute to the 
first integral on the r.h.s. of eq. ( U.l|) , because that integral is of the same isospin 



structure as the one occurring in the Olsson sum rule, but it does generate a 
small negative term via the second integral (—0.02 ±0.01). The net result for the 
asymptotic contributions, 

S as = 1.10 ±0.19 , 

leads to S D +S F +S as = 1.77±0.19. Within the errors, the outcome agrees with the 
numerical value S P = 1.93±0.08 obtained for the l.h.s. of the sum rule (|C.3|) . Note 
that more than half of the r.h.s. stems from the asymptotic region. We conclude 
that our asymptotic representation is valid within the estimated uncertainties, 
also for this sum rule, which converges more slowly than the moments considered 
in appendix B. Since the Olsson sum rule belongs to the same convergence class 
as the one above, we feel confident that our error estimates apply also in that 
case. 



C.2 Asymptotic behaviour of the Roy integrals 

If the imaginary parts of the partial waves with i > 1 are discarded, the Roy 
equations become a closed system for the S- and P-waves. The explicit expres- 
sions for the kernels show that the dispersion integrals over the imaginary parts 
of these waves grow linearly with s, like the subtraction polynomials. Except for 
the contributions from the higher partial waves, the r.h.s. of the Roy equations 
for the S- and P-waves thus grows in proportion to s: 

Ret2-> r , Ret} -> , Ret. 



YIMl ' 1 72M2 ' u 24M2 ' 

„ r , AMI r° , 2 Im tg(s) + 27lmt\(s) - 5 Im t§(s) 



E = 2a° -5al 2L / ds uv ' 1V "\ ^ . (C.4) 

7T i 4 Af| s{s-AMl) v ' 

So, if the coefficient E vanishes, the contribution from the dispersion integrals 
cancels the one from the subtraction polynomial, simultaneously for all three 
partial waves || [82|]. In fact, if the imaginary parts of the higher partial waves are 
dropped and if S is set equal to zero, the Roy equations become identical to those 
proposed by Chew and Mandelstam (see ref. for a detailed discussion). 



The expression for E resembles the Olsson sum rule, where the contributions from 
the S- and P-wave read 

2 AMI f°° 1 2 Im tg(s) + 9 Imt}(s)- 5 Imtg(s) 



2 a u - 5 ai = ^ / ds uw ' ^ ^ + . . . 

7T Jam* s(s-AM^) 
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If only the S'-waves are retained, the Olsson sum rule does imply that E vanishes - 
evidently, this sum rule is closely related to the observation that the linearly rising 
contribution from the subtraction terms must cancel the one from the dispersion 
integrals (section [10]). As is well-known, however, the coefficient of the P-wave 
term in £ differs from the one in the Olsson sum rule. The implications of this 
discrepancy for the Chew-Mandelstam framework are discussed in the references 
quoted above. The family of sum rules derived in appendix |C.1| points in the 
same direction: The imaginary part of the P-wave is tied together with those of 



the higher partial waves - setting these equal to zero leads to inconsistencies [p6 



For the above asymptotic formulae to apply at E ~ 1 GeV, two conditions 
would have to be met: (a) the contributions from the higher partial waves can 
be ignored at these energies and (b) the integrals over the imaginary parts of the 
S- and P-waves are dominated by the contributions from low energies. Unfortu- 
nately, neither of the two conditions is met. The solutions show a pronounced 
structure in the region above the matching point - evidently, we are not dealing 
with the asymptotic behaviour there. The numerical value of £ is negative: The 
integral in eq. ( |C.4| ) over-compensates the term 2a° — 5oq. We may lay the blame 
upon the contributions above the matching point - if the integral were cut off 
there, £ would approximately vanish. 

The situation is quite different for the Olsson sum rule, which does not rely on 
low energy approximations but represents the exact version of the condition that 
must be obeyed by the two subtraction constants for the scattering amplitude 
to have the proper asymptotic behaviour. In that case, the coefficient of the 
P-wave is three times smaller - the region above the matching point plays an 
essential role in bringing the sum rule into balance. The numerical evaluation in 
section [11] shows that even those from the region above 2 GeV are significant. 
The rapid growth of the driving terms indicates that the higher partial waves 
become increasingly important as the energy rises - it is clear that the asymptotic 
behaviour of the partial wave amplitudes cannot be studied on the basis of the 
S- and P-wave contributions to the r.h.s of the Roy equations. 

We conclude that, at the quantitative level, the above simple mechanism 
cannot explain why, for suitable values of and Oq, our solutions remain within 
the bounds set by unitarity. For an analysis of the behaviour above the matching 
point that neither discards the higher partial waves, nor relies on low energy 



dominance, we refer to sections 10 and 11 



D Explicit numerical solutions 

In this appendix, we make available our explicit numerical solutions of the Roy 
integral equations. We proceed as follows. For a few tens of pairs (a{j, ofy in the 
universal band (see fig. |7]), we have constructed the three lowest partial waves at 
2M n < < 0.8 GeV. As we explain in the main text, we parametrize the phase 
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Table 7: Polynomial coefficients for Roy solutions. 
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shifts Sj of the solutions in the manner proposed by Schenk |65 . 

tan^ = ^/l-^S ^ {A\ + Bk 2 + Ctf + D^} (^^f) , (D.l) 

Each solution of the Roy equations corresponds to a specific value of the 3x5 
coefficients in this expansion, 

We approximate these by a polynomial of third degree in the scattering lengths 
a[] and a 2 ,. In terms of the variables 



a a 2 
« = - 1 , v = - 1 , p = 0.225 , p 2 = -0.03706 , 

P0 P2 



the numerical representation for the coefficient B®, for instance, reads 

±? = 2^ + £ 2 M + ,23 V + Z4 U + Z5V + Z e UV + Z 7 U + Z$U V + ZgUV +-210^ . 

The 15 x 10 numbers Zi, . . . ,z w for the coefficients Aq, B®, . . . , s% are displayed 
in table [7], in units of M 2 . 

E Lovelace- Shapiro- Veneziano model 

In this appendix, we describe the model used to illustrate the basic properties 
of the Regge poles occurring in the asymptotic representation of the scattering 
amplitude [^, [84], |85| . In this model, the 7T7T scattering amplitude is taken to be 
of the form 

A(s, t, u) v = Ai $(a s , a t ) + Ai $(a s , a u ) + A 2 $(a t , a u ) , 
where 0) is closely related to the Beta-function, 

1 ' PJ r(i -«-/?) ' 

and a s represents a linear Regge trajectory, 

a s = «o + ckis . 

At fixed t, the function $(a s , at) shows Regge behaviour when s tends to infinity: 

$(a„at) -> (-a s ) Q T(l - a t ) . 
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At the same time, the representation (1 — at > 0) 

$(a s , a t ) = (1 - a s - a t ) B(l - a s , 1 - a t ) (E.l) 

/ 1 ~ a t (a t + 1) •••(<** + ra- 1)1 
= (1 " as " at) l^ + S n!(n + l-a.) j 

shows that the amplitude may be expressed as a sum of narrow resonance con- 
tributions, with mass 

Ml = (ai) _1 (n-ao) , n = 1,2,... 

The coupling constants Ai,A2 may be chosen such that the amplitude does 
not contain resonances with 1 = 2. For this condition to be satisfied, the corre- 
sponding s-channel isospin component 

T 2 (s, t) v = 2 Ai $(a t , au) + (Ai + A 2 ) ($(a s , a t ) + $(a s , a u )) 

should be free of poles in the physical region of the s-channel. This requires 

A2 = ~ Ai = \\ , 

so that the amplitude takes the form 

A(s, t, u) v = -\\ a t ) + $(a S) a u ) - $(a t , , 

T (s,t)y = -|A{3$(« s ,« t ) +3$(a s ,« u ) -$(a t ,a u )} , (E.2) 
T 1 (s,t) y = -A{$(a s ,ai) - §(a s ,a u )} , 
T 2 (s,t) y = -A$(a 4 ,a u ) . 

In the chiral limit, where the Mandelstam triangle shrinks to the point s = 
t = u = 0, the amplitude must contain an Adler zero there. Indeed, the factor 
1 — a s — a t generates such a zero if a = |. Hence the deviation of a from | must 
be of order M 2 , so that a s — | represents a quantity of order p 2 . At leading order 
of the low energy expansion, the behaviour of the amplitude therefore represents 
the first term in the expansion around the point a s = a t = a u = ^, which in 
view of r(|) = ^pa yields 

A(s,t,u) v = n\(a s -l) + 0{p 4 ), (E.3) 

This does have the structure of the Weinberg formula, provided the intercept «o 
is chosen such that a s passes through the value | at s = M 2 , i.e. 



a = \ - ai Ml . 

The lowest levels of spin 1, 2, 3, 4 indeed occur on an approximately linear 
trajectory with this intercept: Fixing the value of the slope with M p , 

a, = J (M 2 - Ml)' 1 , 
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the model predicts 

M h = 1319 (1275) MeV, M P3 = 1699 (1691) MeV, M h = 2008 (2044) MeV, 
where the numbers in brackets are those in the data tables [fTOfl. 



The representation (E.l) shows that for s > 4M^, t < 0, the imaginary part 
of &(a s ,a t ) consists of a sequence of 5-functions: 



Im$(a s ,a t ) = -ir ^ R n {oc t )5(a s - 



n 



n=l 



T(a t + n) 1 

w = r#w = F^ ai(ai+ ) ''' (ai + ri ' • 

For the imaginary part of the s-channel isospin components, we thus obtain 

o \ oo 

Im T°(s, t) v = + R n(ctu)} S(s - M 2 n ) , 

2 «i „=i 

A 7T °° 

ImT 1 ^, t) v = — Y,{M<Xt) - Rn(a u )} 5(s - Ml) , 
lmT 2 {s,t) v = 0, 

with w = AMI - M l- 

We may then read off the imaginary parts of the partial wave amplitudes by 
decomposing the polynomial R n (a) into a Legendre series:^ 

n ( It \ 

o \ oo 

Imt° e (s) v = -^{l + (-lY}j:rneS(s-Ml) , 
64 ax ^ 

\ 00 
Im4(«)v = - (-1)'} E ^(s - Ml) , 

lmt 2 e {s)v = . 
On the leading trajectory, the coefficients are 

r nn = — — (Ml - 4 M^) n . 

Comparison with the narrow width formula ( |B.13| ) shows that the model 
predicts the width of the various levels as^ 

^= ^f" n U Ml-4Ml)K (E.4) 
62 tc «i M„ 



14 In the case of tg (s) , the sum over n only starts at n = 1. 

15 The formula reproduces the numerical results in Table I of ref. |85| , if the parameter values 
are adapted accordingly (a = 0.48, a x = 0.9 GeV" 2 , T p = 112 MeV, M p = 764MeV). 
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where u 1 depends on the isospin of the particle: u° = 3, u 1 = 2, u 2 = 0. In 
particular, the result for the width of the p reads 

r '=967M^- 4M - )f - < E ' 5 > 
Fixing the coupling constant with the experimental value T p = 151.2 MeV, we 



obtain X/32n = 0.728. The formula ( |E.4|) then predicts 

= 130 (157) MeV, r™ = 51 (51) MeV, = 46 (35) MeV, 

where the numbers in brackets are again taken from the data tables [ff(|. This 
shows that the model does yield a decent picture, not only for the masses but 
also for the widths of the particles on the leading trajectory. 

In addition to the levels on the leading trajectory, the model, however, also 
contains plenty of daughters, with a rather strong coupling to the 7T7r-channel. 
For the states on the first daughter trajectory, for instance, equation ( |F.4[ ) yields 
= 783 MeV, = 154 MeV, r™ = 113 MeV, = 42 MeV, etc. The scalar 
daughter of the p is particularly fat. 

It is clear that an amplitude that describes all of the levels as narrow reso- 
nances fails here. Unitarity implies the bound 

/ dslmt° (s) Jl -AMl/s < M 2 -AMI. 

This condition is violated for M < 1.3 GeV. Also, if the intercept of the leading 
trajectory is fixed with the Adler condition as above, the scalar daughter of the f2 
is a ghost: The formula ( |E.4|) yields a negative decay width |85|. In this respect, 
the model is deficient - as witnessed by the life of even royal families, the decency 
of a mother does not ensure that her daughters behave. 

The problem also shows up in the S-wave scattering lengths: Chiral symmetry 
relates the coefficient of the leading term in the low energy expansion ( |E.3| ) to 
the pion decay constant, 

7rA«i = — . (E.6) 

If the coupling constant A is fixed such that the model yields the proper width for 
the p, the l.h.s. of this relation exceeds the r.h.s. by a factor of 1.7. Accordingly, 
the prediction of the model for a[j exceeds the current algebra result by about this 
factor. In the vicinity of threshold, the behaviour of the amplitude is determined 
by the properties of the function (fi(a,j3) for a ~ /3 ~ ^. There, the first term in 
the series ( |E.1| ) accounts for about two thirds of the sum. The spin 1 component 
of this term is due to p -exchange, while the spin part arises from the scalar 
daughter of the p. By construction, the former does have the proper magnitude. 
The S-wave scattering lengths are too large because the scalar daughter of the p 
is too fat. 
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As was noted from the start f85 |, the model is not unique. To arrive at a more 
realistic model, we could add extra terms that domesticate the daughters and 
leave the leading trajectory and the position of the Adler zero untouched. Note, 
however, that the number of states occurring in the Veneziano model corresponds 
to the degrees of freedom of a string, while the spectrum of bound states in QCD 
is the one of a local field theory, where the number of independent states grows 
much less rapidly with the mass. Modifications of the type just mentioned can at 
best provide a partial cure. In particular, these do not remove the main deficiency 
of the model, lack of unitarity. 
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